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1 Introduction 



La dynamique des populations est au cœur de l'interface entre systèmes 
dynamiques et biologie. Ainsi, l'un des modèles biologiques les plus simples 
et les plus importants — le modèle logistique — correspond à la dynamique 
des polynômes quadratiques, dont l'étude mathématique est des plus inté- 
ressantes. Ce seul cas montre combien il est aisé d'obtenir un comportement 
complexe sans introduire beaucoup de complexité dans le modèle. 

Mais le modèle logistique est bien trop peu réaliste pour que sa complexité 
dynamique puisse être directement interprétée dans le cadre d'une popula- 
tion biologique réellement observée. D'un point de vue mathématique, les 
systèmes de dimension 1 présentent un nombre limité de dynamiques pos- 
sibles. Il est donc intéressant de considérer des modèles de dimension supé- 
rieure tels que l'application de Hénon (qui est une petite perturbation de la 
dynamique d'un polynôme quadratique, en dimension 2), et qui sont encore 
mal compris du point de vue théorique. 

Le pas suivant dans cette démarche consiste en l'étude de systèmes dy- 
namiques de dimension infinie, qui seront de «petites perturbations» des 
modèles précédents, i.e. le modèle logistique. En effet, si l'on veut intégrer le 
phénomène biologique de maturation des jeunes individus, il est nécessaire 
de considérer la fonction d'évolution de l'effectif en temps continu, et non 
seulement sa valeur à un instant donné, ce qui donne un système dynamique 
de dimension infinie (ou de grande dimension, si l'on discrétise ce système). 
Un autre phénomène intéressant à considérer est l'infiuence des rythmes sai- 
sonniers sur un tel système, lorsqu'il se combine avec cet effet de retard induit 
par le temps de maturation des jeunes. Le modèle que nous considérons com- 
bine ces deux effets avec une forme de densité-dépendance un peu différente 
de celle du modèle logistique. 

Nous commencerons par définir le modèle étudié, tel qu'il a été énoncé 
dans |YB98] . puis sous une forme légèrement modifiée, en motivant celle- 
ci aussi bien par des raisons biologiques que des raisons de simplicité pra- 
tique. Nous verrons ensuite ce que l'on peut montrer simplement par une 
étude théorique a priori, point de départ d'une étude future plus approfon- 
die (mais surtout bien plus difficile). La dernière et plus importante partie 
de notre étude sera consacrée à l'analyse des résultats de simulations numé- 
riques, en vue de comprendre l'infiuence des paramètres sur la dynamique 
du système et d'analyser plus finement un des attracteurs étranges que nous 
avons pu observer. Ce travail est bien sûr loin d'être complet, et se veut sur- 
tout être une introduction (et une motivation) pour de futurs travaux, aussi 
bien mathématiques que biologiques. 
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2 Description du modèle 

2.1 Modèle initial continu 

Le modèle suivant, défini dans |YB98j . décrit l'évolution temporelle d'une 
population de campagnols. 

i-Ai 

N{t)= S{a)mp{t- a)N{t- a)m{N{t- a))da (2.1) 

- t est le temps (en années), 

- la population active (z.e. d'âge supérieur à Aq), 

- Aq l'âge de maturation, 

- Al l'âge maximal, 

- 5* le taux de survie, 

- rrip le paramètre de saison (décrit une probabilité de reproduction en 
fonction de la saison), 

- m{N) le taux de fécondité individuel annuel pour une population active 
de taille N. 

On a choisi des formes simples pour les fonctions S, rUp et m : 

S{a) = 1 - ^ (2.2) 

fo si < t < p mod. 1 , , 

m,{t) = \ . - -/ (2.3) 
I 1 SI p < t < 1 

, , 1 nin si < 1 , , 

m{N) = \ \ - (2.4) 
\mQi\ ' smon 

Cependant, pour éviter des artéfacts dûs à la non-régularité du système, 
il nous a semblé préférable de remplacer m et rUp par des fonctions un peu 
plus régulières. 



2.2 Lissage de la fécondité 

La fonction m définie par (12. 4p est continue mais pas C^. Il y a une forte 
rupture de pente à la valeur critique N = 1. On pourrait aisément «recoller» 
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les deux parties de la courbe pour obtenir une fonction C°°, mais cela ne 
serait pas très pratique pour les simulations numériques. Nous utilisons ainsi 
une parabole intermédiaire qui rend m C^. 



niosi N < Ni^^ 
m{N) = {mox {A^ + B^N + C^N'^) si A^i,^ < < A^a^ 
moN-^ si A^2,7 < N 



(2.5) 



Pour garder un modèle réaliste et suffisamment proche du modèle initial, 
il faut conserver la décroissance de la fécondité h-> m{N), et s'assurer que 
la parabole rejoint les valeurs extrêmes rriQ et îïIqN^"' suffisamment près de 

= 1. Pour définir complètement les paramètres A^^i ,^, A^, B^, et A"2,^, 
on impose également m{N2^^) = mo/2. Cette dernière contrainte permet 
de limiter la zone intermédiaire, ce qui simplifiera notamment les calculs 
explicites de la section [31 On a ainsi les conditions suivantes : 



B^ + 2C^Ni^^ 

2,7 



B^N2^^ 



1 



1 

2 
1 



-7-1 



2,7 



7 



2A^. 



2,7 



(2.6) 



De plus, on veut < pour garantir la décroissance de m. On choisit 
donc : 



'^2,7 


= 2^/^ 




C'y 






8 X 41/7 




^ ^7 




4 


B^ 


= 2-1/^ X 


4 


^1,7 


= 2^/^ X (^1 


2' 

7y 



(2.7) 



La figure [T] représente les deux fonctions m — continue et — pour 
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8,25. 
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Comparaison des fonctions de fécondité ; = 50 ; y = 8.25 



- - fécondité continue 

— fécondité 



N = effectif mature 



Figure 1 - Comparaison des fonctions de fécondité. 



2.3 Lissage du facteur saisonnier 

La fonction ntp définie par (12. 3p n'est pas continue, il est légitime de 
vouloir considérer un facteur saisonnier un peu plus régulier. On a choisi, 
arbitrairement, de le rendre en effectuant le passage de à 1 à l'aide d'un 
cosinus. Pour cela, on ajoute un paramètre e qui est la durée du printemps 
et celle de l'automne. La durée de l'été est désormais 1 — p — e et non plus 
1 — p. Prendre e = ramène bien sûr au cas précédent. 



ri(l + cos(7rx(i + i))) 



rUp^^it mod. 1^ 



\ (1 + cos (tt X 

1 

i(l + cos(7r X 



-m 
m 



si < t < e/2 
si e/2 < t < p - e/2 
si p - e/2 < t < p + e/2 
si p + e/2 < t < 1 - e/2 
si 1 - e/2 < t < 1 

(2.8) 

Cette définition n'étant valable que lorsque e < p < 1 — e, on posera 
e = min(p, 1 — p) lorsque ce n'est pas le cas a priori. 

La figure |2]représente mp^^ pour p = 0,41 et deux valeurs de e. Remarquons 
enfin que l'on aurait également pu rendre C°° dès que e > en utilisant 
autre chose qu'un cosinus. Le choix que nous avons fait tient compte de la 
simplicité des calculs numériques futurs. 

Pour les détails concernant les simulations numériques, voir le paragraphe 
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Régularisation du facteur saisonnier (p = 0.41) 




Saisons non-continues : e^^^ = 
Saisons :e , =0.1 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

t en années 



Figure 2 - Comparaison des fonctions rrip^^. 



EU 



en annexe. 



2.4 Aspects biologiques 
2.4.1 Études antérieures 

Espèces concernées Le modèle que nous venons de décrire a été élaboré 
en relation avec l'étude de la dynamique de certaines populations de petits 
rongeurs. Ceux-ci se caractérisent en effet par un fort investissement dans 
la reproduction (portées importantes et maturité sexuelle très rapide) et de 
grandes variations annuelles de la taille de la population. Plusieurs espèces 
de campagnols ont ainsi été étudiées, notamment le campagnol roussâtre 
Clethnonomys glaerolus (figure [3]; (YM lCVS+n2[ lYSHP.Tni] ) et Microtus 
townsendii [LYOlj . 

C'est plus particulièrement une population de Microtus epiroticus^ in- 
troduite accidentellement il y a une cinquantaine d'années dans l'archipel 
arctique de Svalbard (en Finlande, dans le Spitzberg), qui est visée par ce 
modèle. Ces campagnols possèdent en effet une fécondité est extrêmement 
élevée pour des mammifères. De plus, cette espèce étant menacée d'extinc- 
tion, l'étude de sa démographie permettrait également de mieux la protéger. 
La figure m montre ainsi de grandes fluctuations de population, et des effectifs 
minimaux très faibles, de l'ordre de quelques individus. 



1. sibling vole en anglais. 



5 
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Figure 3 - Campagnol roussâtre (bank vole, Clethrionomys glareolus). Ce 
petit rongeur, de la famille des Microtinés, vit dans les forêts tempérées et 
se nourrit principalement de graines. 
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Figure 4 - Microtus epiroticus à Svalbard : nombre d'individus capturés 
entre 1991 et 2002. (Source : NINA and University of Troms0. 2003 © Sta- 
tistics Norway). 



Mécanismes envisagés Plusieurs causes possibles à ces phénomènes ont 
été étudiées. Il a été montré expérimentalement qu'une augmentation de 
la quantité de nourriture disponible augmente la densité mais n'a pas d'in- 
fluence sur la densité-dépendance (C. glaerolus, Finlande |YSHPJÔT] ). Plus 
précisément, la nourriture disponible est liée aux variations inter-annuelles, en 
jouant sur la survie de l'année suivante ; on observe des variations saisonnières 
de taille comparable aux variations inter-annuelles et qui se caractérisent par 
une survie plus faible au printemps (C. glaerolus, Belgique |CVS"'"02] ). 

Le facteur climatique semble lui aussi jouer un grand rôle : une compa- 
raison entre M. epiroticus à Svalbard et Chionomys nivali^ dans les Alpes 
françaises indique une relation entre un environnement stable et un faible tur- 
nover I [Yl99]. En effet, la population vivant dans l'Arctique, où les conditions 
hivernales sont très variables, a une reproduction extrêmement rapide, tandis 
que la population alpine, dont l'environnement est stable, se reproduit peu 
et a une forte survie. On a ainsi essayé d'inclure une stochasticité environne- 
mentale et démographique dans les modèles afin d'expliquer les fluctuations 
de population observées (C. glaerolus, Alpes |YSOO] ). 

Un phénomène important pourrait être relié à ces facteurs environnemen- 
taux : la plasticité de l'âge à maturité. Il a ainsi été montré que les femelles 
nées en début de saison de reproduction ont avantage à se reproduire ra- 

2. snow vole en anglais. 

3. i.e. forte survie et faible fécondité 
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pidement, malgré le coût élevé d'une reproduction précoce (M. townsendii, 
Vancouver |LY01] ) : les femelles naissant plus tôt peuvent se reproduire avant 
le fin de la saison de reproduction, augmentant ainsi la contribution de leur 
mère à la population totale. 



2.4.2 Paramètres du modèle 

Les paramètres de croissance et de reproduction des deux populations^ 
ont été évalués dans [ YIS93J . A partir de ces conclusions, nous pouvons jus- 
tifier le modèle et le choix de paramètre «typiques». 



Densité-dépendance La forme particulière de densité-dépendance se jus- 
tifie car le facteur limitant est ici le nombre très restreint de sites de repro- 
duction. Ainsi, seule la population (femelleEl) mature doit entrer en ligne de 
compte, et en cas de surpopulation, les quelques sites disponibles perdent 
beaucoup en qualité. Ceci est confirmé par les conclusions de |YSHPJÔT] . qui 
indiquent dans le cas de C. glaerolus une densité-dépendance plus forte en 
été, à cause de la maturation des femelles. Nous choisissons de ne pas faire 
dépendre 7 de la saison puisque la variation observée est liée à la prise en 
compte de la densité totale, et non de la densité de femelles matures comme 
nous le faisons ici. 



Saison de reproduction La période de reproduction correspond à la sai- 
son de croissance des plantes, c'est-à-dire du printemps à l'automne. Nous 
considérons avec ce modèle un climat parfaitement prévisible, identique d'une 
année sur l'autre. Il s'agit de savoir si l'on peut observer un comportement 
chaotique dans des conditions parfaitement stables. 

Sa durée varie donc selon les lieux. À Svalbard comme dans les Alpes, 
celle-ci dure de 3 à 4 mois (Juin à Septembre) |YI99] . i.e. p ~ 0,7. Dans 
la plupart des environnements tempérés, celle-ci est beaucoup plus longue. 
Ainsi, en Belgique, C. glaerolus se reproduit en général entre la deuxième 
semaine d'Avril et la fin du mois d'Octobre |CVS+n2l , mais peut varier de 
telle sorte que l'on a 0,35 < p < 0,45. 



Fécondité Un élément important du modèle est la forte fécondité des cam- 
pagnols. Pour ajuster les paramètres du modèle, nous avons besoin d'évaluer 

4. M. epiroticus à Svalbard et M. arvalis en Finlande. 

5. Comme souvent en dynamique des populations, seule la population femelle est consi- 
dérée dans la mesure où elle est le facteur limitant de la reproduction. Il n'est intéressant 
de considérer les mâles que si la femelle a peu de chances de rencontrer un mâle (par 
exemple si la densité de population est faible), ou bien si le sex-ratio est loin de 1 : 1. 
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le nombre de jeunes femelles par femelle mature et par an, en l'absence de 
densité-dépendance. 

Pour M. epiroticus |YIS93] . la période de gestation est de 20 jours, et 
la taille des portées varie de 4,5 à 7, celle-ci augmentant pour une même 
femelle au fur et à mesure de ses reproductions. Le sex ratio est très proche 
de 1 : 1. Dans la mesure où une femelle peut se reproduire pendant la période 
d'allaitement, on en déduit une valeur maximale itiq 18 x 6 x 0,5 = 54. 

Dans le cas de M. townsendii [LYOl] , on observe de 5 à 6 portées par 
saison de reproduction (celle-ci durant de Mars à Novembre, soit environ 8 
mois : p ^ 0,3), chacune comportant 5 à 8 individus. On a donc dans ce cas 
mo/3 ~ 5 X 6,5 x 0,5 soit ttiq ~ 24. Il n'est pas étonnant de constater que 
cette valeur est bien inférieure à celle de M. epiroticus, qui représente un cas 
extrême parmi les mammifères. 

On peut donc prendre comme valeur rrio = 50, mais une fécondité légère- 
ment inférieure serait sans doute plus réaliste. 

Âge de première reproduction L'âge de maturité des femelles est su- 
périeur à 17 jours, auquel il faut rajouter la durée de gestation, soit 20 jours 
supplémentaires |YIS93| . L'âge de première reproduction de M. epiroticus est 
donc au minimum 37 jours, i.e. Aq > 0,10. En moyenne, on observe plutôt 
une première portée â un âge d'environ 50 jours, i.e. Aq ^ 0,14. 

Une telle précocité ne se retrouve pas chez les autres Microtinés, â l'excep- 
tion de M. arvalis en Finlande. On utilisera donc 0,1 comme valeur minimale, 
tandis que Aq = 0,20 (correspondant à 75 jours) est plus classique. 

Survie Pour M. epiroticus â Svalbard, la survie hivernale est de l'ordre de 
0,1 (et très variable selon les années), et la survie estivale 0,85 |YI99| . Le 
taux de mortalité par année est donc de 0,046 en hiver et 0,52 en été. L'âge 
maximal Ai est toujours inférieur â deux ans. 

Chez C. glaerolus, en Belgique, les variations saisonnières de la survie 
ont été étudiées plus en détails |CVS"'"02] . Le taux de survie est plus élevé 
en hiver (0,95 par semaine, soit 0,07 par an) qu'en été (0,90 par semaine, 
soit 0,004 par an). Elle est également légèrement différente chez les femelles 
matures et immatures. 

Le taux de mortalité est supposé dans le modèle constant avec l'âge, et 
indépendant des saisons, ce qui est loin d'être le cas en général. L'âge maximal 
de 2 ans est également une légère sur-estimation de ce qu'il est en réalité. 
Pour plus de réalisme, la fonction de survie est sans doute l'un des premiers 
éléments du modèle â modifier. 
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2.4.3 Problêmes posés 

L'étude de ce modèle n'a pas pour but de faire des prévisions précises 
concernant l'avenir de la population de Microtus epiroticus à Svalbard. Nous 
nous efforcerons de considérer des paramètres réalistes pour de petits ron- 
geurs, pas nécessairement M. epiroticus. 11 s'agit surtout d'effectuer un travail 
théorique dans un cadre assez simple, afin de savoir si l'on peut observer une 
grande variabilité d'effectifs (voire une dynamique chaotique) dans un mo- 
dèle complètement déterministe, dans un environnement régulier. Et si oui, 
quelles sont les facteurs biologiques déterminants (fécondité, âge de maturité, 
durée de l'hiver) ? Un autre objectif est de bien comprendre ce modèle très 
simple avant de le complexifier en introduisant d'autres mécanismes pouvant 
jouer un rôle dans la dynamique de cette population, parmi ceux que nous 
avons évoqués précédemment. 



3 Étude théorique 

Les résultats de cette partie proviennent principalement de |YB98] . où les 
fonctions mp et m considérées étaient données par (12. 3 p et (12.4p . Nous avons 
considéré ici un cadre un peu plus général, valable pour les fonctions irip et 
m «lissées». 



3.1 Définition du système dynamique 

Pour to G M/Z, notons Yt^ l'ensemble des fonctions continues sur 
[—Al] 0] à valeurs positives et vérifiant la condition 

/•Al 

iV(0) = / S{a)N{-a)m{N{-a))mp{tQ-a)da. 

Le système dynamique défini par (12. ip est donné par le semi-groupe {Ts)s>o '■ 

T'{t, N) = {t + s(mod. 1), N^) (3.1) 

N{s - a) si < s < a < Al 
* /a ^ S{b)N{s — a — b)'m{N{s — a — h))mp{t + s — a — h)dh sinon 

(3.2) 

Ceci est bien défini pour < s < Aq, on l'étend à s > par la propriété de 
semi-groupe. L'espace des phases est alors 

ytt = {{t,N)/t e M/Z, G Yt}. 



Nfi-a) 
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Dans la suite, on écrira parfois A^* au lieu de Nf lorsque cela ne crée pas de 
confusion. 

On s'intéresse à l'application : Yq ^ Yq donnant l'évolution de la 
population d'une année sur l'autreH A N E Ytg donnée, on associe ainsi 
une unique fonction continue définie sur [— y4i;+oo[ prolongeant et 
compatible avec T (c'est-à-dire la solution de l'équation f l2.ip ). 

On munit Y^ de sa topologie naturelle, induite par la topologie produit 
sur M/ZxC([— Ai; 0]), l'espace C{[—Ai] 0]) des fonctions continues de [—Ai, 0] 
dans R étant muni de la topologie de la convergence uniforme. Il découle alors 
de la continuité uniforme de — > x m{N) la propriété suivante : 

Proposition 3.1. Pour tout s >0, est un opérateur continu Y^ YK 

De plus, N ^ N X m{N) étant f^j-Lipschitzienne (la constante Kf peut 
être calculée explicitement en fonction des paramètres du modèle) , l'applica- 
tion T"^ : Yt ^ Yt+s est i^'-Lipschitzienne, avec K = max(l, [Ai — Aq) x Kj). 
Ceci découle directement de la définition de N^. La constante K dépend 
uniquement des paramètres du modèle, et pas de i G R/Z. 

3.2 Existence d'un attracteur 

Nous allons montrer que pour des valeurs raisonnables des paramètres, 
un tel système dynamique possède un attracteur, ce qui nécessite plusieurs 
lemmes techniques. Nous n'utiliserons pas les formes explicites des fonctions 
rrip et m (pour rester généraux, dans la mesure où celles-ci pourraient être 
modifiées ultérieurement), mais uniquement les hypothèses suivantes : 

si N <1 (3.3) 

si N>1 (3.4) 
Vt (3.5) 

mp{t) = 1 sur un intervalle de longueur 1 — p — e. (3-6) 

On voit aisément que les fonctions définies par (12. 8p et (12. 5p (resp. (12. 3p 
et (12. 4p ) vérifient ces hypothèses. Le symbole A est employé ici et dans la 
suite à la place de min, de même que V signifie max. 

6. Le choix de to = est arbitraire, on pourrait aussi bien considérer : Ytg — Yt^, 
qui donnerait le même type de dynamique. 



niQ > m[i\) > 



moN-^ > m{N) > y-^ ^ > 
1 > m„(t) > 
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Posons 

f-Ao+l 

Cq : = S{a)da (3.7) 

J Ao+p+e 

= {l-p-e)[l — J. (3.8) 

Nous nous plaçons désormais dans le cas où les paramètre vérifient les condi- 
tions suivantes : 

7 > 1 (3.9) 
Al > {2Ao) V (Ao + 1) (3.10) 
co X mo > 2 (3.11) 
p + e<l (3.12) 

Remarque 3.1. Ces conditions sont très raisonnables, et toujours vérifiées au 
cours des simulations que nous avons faites. En effet, si on impose p + e < ^, 
Al > 2, Aq < 1/2, 7 > 1, alors cq > et donc mo > 15 suffit pour 
satisfaire fl3.1ip . Il n'y a donc pas à s'inquiéter du manque de finesse de cette 
majoration. 

Lemme 3.2. Soit to G M/Z, G Yt^^. On a alors, pour tout < s < Aq : 

N{s) < iV^ax := mo Y (l - ^ 
Démonstration. On a toujours Nm{N) < niQ, d'après (13.3p . (13.4p . et car 



7 > 1. Comme de plus rup < 1 (13.5p . on a 



N{s) < mo / S{a)da = 

'Ao 



□ 



Lemme 3.3. Soit N G Yf^^ telle que N < A^max- On a i{N) = min[_yi^.o] > 
0. 

1. Sîi{N)< N}^l2, alors N{s) > ^i{N) pour < s < Ao. 

2. Si i{N) > N^-J, alors N{s) > ^N^^J Pour < s < Aq. 

Démonstration. Commençons par montrer que A^m(A^) > ^x[i[N) A N^J) 
sur [-Ai;0] : 



12 



si < 1, 



et si > 1, 



Nm{N) >Nx^ 

> .(AT) X —, 



Nm{N) > — X N^-^' 

2, (3.14) 
> !^ X iv^-^. 

— 2 max 



Par conséquent, 



N{s) > ^ X (i(Ar) A Ar^;;j) x / ' S{a)mp{to + s - a)da. 

En fonction de la valeur de to + s, on peut trouver un sous-intervalle de 
[^o; ^0 + 1] C [Aq] Al], de longueur 1 — p — e sur lequel rripito + s — ■) vaut 1. 
La dernière intégrale est donc minorée par la même intégrale restreinte à ce 
sous-intervalle, qui est plus grande que cq car S est décroissante. Ceci achève 
la preuve du lemme. □ 

Corollaire 3.4. S'z 7 > 1, ^^^^ > 1, N & Yt^, alors pour s assez grand 
(dépendant de N), on a : 

^-^N'^i: < N'ia) < Nrna., Va G [-Ai;0]. 

Remarquons que l'on peut remplacer la condition fl3.1ip par CorriQ > 1 
dans le cas où m est définie par (12. 4p (cf. |YB98] ). La constante 2 a été choisie 
arbitrairement dans l'opération de lissage de m, celle-ci pourrait être prise 
plus proche de 1 sans difficulté supplémentaire, mais toujours strictement 
supérieure à 1. 

Lemme 3.5. Soit N G yi,,. Posons L = mo — Alors, si < sq < 
Si < Aq, on a 

\N{si) -N{so)\ <L\si-so\. 

Démonstration. Remarquons tout d'abord que nip < 1, < S < 1 et 

|5'(si - u) - S{so -u)\< A];^ \si - Soi ■ 

Plaçons-nous dans le cas où Si — Ai < sq — Aq (c'est vrai car on a supposé 
y4i > 2 X Aq), et écrivons la définition de N[si) en fonction de N. 
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|iV(si)-iV(so)| 



+ 
+ 



si-Ai 

S{sq — u)N{u)m{N{u))mp(to + u)du 

sa- Al 

[S{si — u) — S{sq — u)]N{u)m{N{u))mp{tQ + u)du 

si-Ao 

S{si — u)N{u)m{N{u))mp(to + u)du 

so-Aq 



< mo |si - Sol + mo 1 1 - ^ 1 |so - Si| + mo |si - So| 



< L |so — si| 

Nous pouvons maintenant définir 

Kt, =[n e Yt,- Vs G [-Ai; 0], ComoN^i: < N{s) < N, 
Wso, s, e [-A^■, 0], \N{so) - N{s,)\ <L\so- s. 



(3.15) 

(3.16) 
(3.17) 

□ 



(3.18) 



qui est une partie compacte de Yt^ pour la topologie de la convergence uni- 
forme, d'après le théorème d'Ascoli. 

Les lemmes que nous venons de démontrer peuvent se formuler de la façon 
suivante : 

Proposition 3.6. Soit N G Yq, (T^)s>o ^6 semi-groupe défini par l'équa- 
tion (13. ip . On se place dans les conditions précédemment énoncées pour les 
différents paramètres du modèle. 

1. Si N & /Co, alors N'"^ G JCg pour tout s > 0. En particulier {Kq) C /Cq- 

2. En général, il existe Sq > (dépendant de N) tel que G ICg pour 

tout s > Sq. 

Démonstration. 1. N'"^ G Yg par définition, A^^ < N^ax d'après le lemme 
13. 2[ A^^ reste L-lipschitzienne d'après le lemme 13.51 La partie 2 du 
lemme [3731 donne la minoration, en utilisant que A^ G Yq. 

2. On utilise le corollaire 13.41 pour montrer l'existence de Sq, le reste de la 
preuve étant identique. 

□ 
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U attracteur du système dynamique (F^, (T*)s>o) est défini par : 

A = {(t, N)/t e M/Z, N e ^t} (3.19) 
avec = Pi T"(/Ct) (3.20) 

n>0 

La propriété suivante justifie l'appellation d'attracteur pour A. 

Proposition 3.7. 1. A est une partie compacte de YK 

2. Pour tout s > 0, T*(A) = A. 

3. Pour tout voisinage U de A, et toute condition initiale (0, A^), G Yq, 
il existe sq (dépendant de N et U) tel que T'^{0,N) G U pour tout 

s > Sq. 

Démonstration. 1. D'après la continuité de et la compacité de ICt, At 
est compact pour tout t > 0. De plus, > Id uniformément sur 

(c'est une conséquence du lemme [ÏÏ75| car on a alors ||T^(A^) — A^||oo < 
Le). Par conséquent, A est compact. 

2. Par construction, T^{{t} x A^) = {t + s} x At+s pour tous s,t G M/Z, 
d'où T%A) = A pour tout s > 0. 

3. On peut supposer N G JCq d'après la proposition 13.61 La suite T"(0, A^) 
est alors contenue dans )C = {{t, N),t G M/Z, A^ G ICt} qui est compact. 
Tout point d'accumulation de cette suite est nécessairement dans A, et 
donc T"(0, A^) G U pour n entier assez grand. 

De même, en considérant la suite (T""'(0, N))^^^ avec a > réel quel- 
conque, on montre que T""'(0, A^) G U pour n > n{a, U) entier. Munis- 
sons de la distance 

diis,N),{t,N)) = \s-t\ + \\N-N\U 

qui engendre bien la topologie de y" précédemment définie. Puisque les 
éléments de /Cq sont L-lipschitziens, on a pour tous s,t >0, 

d{T'{0,N),T\0,N)) = \s-t\ + ||A^" - A^*||oo (3.21) 

< |s-t| + ||ÎV(s + -)-ÎV(t + -)||oo (3.22) 
<\s-t\{l + L). (3.23) 

De plus, comme A est compact, il existe e > tel que A C A^^^ C U où 
l'on a noté A*^^) l'épaississement de e de A (c'est-à-dire l'ensemble des 
points situés à distance < e d'un point de A)0. 

7. on recouvre A par des boules contenues dans U, un nombre fini suffit par compacité, 
e est alors le min des rayons de ces boules. 
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Prenons a = 2{i+l) ' ^1*^^^ pour tout t > n(a, A^*^/^^) réel, T*{0,N) G 
A*^^) C t/ ce qui achève la preuve. 

□ 

Il est alors aisé de faire le lien avec les définitions données IC.14I et IC.15I 
données en annexe, sous la forme du corollaire suivant. 

Corollaire 3.8. Le compact A défini par l'équation fl3.19p est un attracteur 
pour le système dynamique ((T'')^>q , y**) . Son bassin d'attraction est YK 

Démonstration. La seule difficulté supplémentaire par rapport à la propo- 
sition 13.71 est qu'il faut montrer l'existence d'un voisinage de A revenant 
tout entier dans lui-même en un temps fini A^. Pour l'instant, nous savons 
seulement que toute condition initiale arrive en temps fini dans un voisi- 
nage donné, mais ce temps peut être arbitrairement grand en fonction de la 
condition initiale dès que l'on est hors de A. 

D'après le 2. de la proposition 13. 7[ il revient au même de considérer le 
système dynamique discret {T^,Y^). Nous nous placerons désormais dans ce 
cas. Soit A^^i > Al un entier et a > tel que (1 — a)^^^ > 1. On définit alors 
l'ensemble 

y = ^TN.y^ [{it,N)/t G M/Z, NeY,,N> N^-^ x ^ x (1 - a)}) 

Par continuité de T^^, c'est un ouvert. Par définition de /C^, il contient 
{(t, N)/t G M/Z, N G /CJ, et donc A. 

Les lemmes [321 O et [33] montrent que T^^^{V) C {(t,A^)/A^ G /CJ, 
d'où T^^^{V) C V. 

On a également 

fl T^{V) C fl T"({(t, N)/N G )Ct}) = A 

n>0 n>0 

Comme de plus A C V et T"(A) = A, on a A = n„>o^"(^)- 

Le 3. de la proposition 13.71 montre que si a: = {t,N) G y, il existe un 
temps to à partir duquel T*(x) G A, et donc nécessairement u{x) C A. □ 



3.3 Modèle non-saisonnier 

Le cas rrip = 1 (p = 0) peut être traité plus en détails, au moins pour de 
petites valeurs de 7. 

Il existe une valeur d'équilibre (z.e. une solution constante en temps 
continu) 
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pourvu que la quantité 



2Ai 



(^1 - Ao) 



S{a)da 



soit plus petite que mg. Si elle est plus petite que lorsque m est définie 
par (12. 5p . on peut réécrire 



eq 



niQ- 



(^ 1 - ^o) 
2Ai 



1/7 



C'est toujours le cas pour des valeurs raisonnables des paramètresH. 
Pour déterminer la stabilité de cet équilibre, définissons 



F(A) 



S{a)e-''''da 

Ao 

^11-^1- 



AJ X^Ai 



-A) A 



-AiX 



(3.25) 
(3.26) 



les valeurs propres de la différentielle!^ à l'équilibre sont les solutions de 

2Ai 



F{X) 



{Al - A^y 



:i-7) 



-1 



On peut alors définir des valeurs 7o(Ao) < 7i(y4o) < 
telles que : 



(3.27) 
< 7fc(Ao) < ••• 



- si 7 < 7o(v4o) := 1 + 2A^"'' l-^(^^^o)| ^! l'équilibre est stable. 

- si 7o(Ao) < 7 < 7i(^o) := 1 + ^^^^ 



2Ai 



\F{—iu2)\ ""^j il y a exactement 



deux valeurs propres (complexes conjuguées) instables. 



- si 7fc_i(v4o) < 7 < 7fc(^o) := : 
ment 2k valeurs propres instables!^. 



2Ai 



il y a exacte- 



e.g. 7 > 1, ^1 > 2, Al > 2Ao, mo > 



9. On considère le système sous la forme ^ = f{N), au voisinage de N^q (on peut 



expliciter /, au voisinage de N^q, en considérant T*^ quand e tend vers 0). Les solutions 

'^'^ — D/n^ h. Si X est valeur propre de la différentielle à 



de la forme Neq + h vérifient 



l'équilibre, un vecteur propre associé h est nécessairement sous la forme h{t) = /loC^*. 
En introduisant cette relation dans (I2.ip . on peut ainsi déterminer l'ensemble des valeurs 
propres. 

10. k paires de valeurs propres complexes conjuguées. 
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Eléments de preuve : Les valeurs propres de la différentielle sont stables si et 
seulement si leur partie rélle est négative, puisque le vecteur propre associé 
est de la forme t — )■ exp(At). Comme c-y est un réel (négatif si 7 > 1), ce 
sont des solutions de l'équation '^{F{X)) = 0. Les valeurs des bifurcations 
correspondant à A imaginaire pur, on s'intéresse à l'équation Q{F{—iu)) = 0, 
avec u > 0. 

On montre alors que l'ensemble de ces solutions peut s'écrire Uq < Ui < 
U2 < ■ ■ ■ avec 

r ^{F{-lU2k)) < '^{F{-lU2k+2)) < • ■ ■ < 
\ ^{F{-lU2k+l)) > ^{F{-tU2k+3)) > • ■ ■ > 

En effet, 

'^{F{~iu)) = sin(Ain) + -(1 - 4^) cos(Aom) + ^— sm{Aou) 
w^Ai u Al u^Ai 

3î(F(-m)) = cos(Ain) + -(1 - 4^) sin(v4oM) + cos(Aom) 
u^Ai u Al w^Ai 

et donc l'ensemble des solutions n'a pas de point d'accumulation (en 0, on 
le vérifie par un développement limité de 'is{F{—iu)) ; ailleurs, cela découle 
clairement de l'analycité de la fonction qui n'est pas identiquement nulle), ce 
qui permet d'énumérer les solutions. Il faut vérifier par un calcul direct que 
les inégalités annoncées sont vraies pour les premières valeurs de k. Pour k 
grand, le terme en 1/u est dominant, et donc '^{F{—iu)) s'annule presque 
en même temps que cos{Aou), et en ce point 3Ç(F(— in)) se comporte comme 
^(1 — ^)sm{AQu). On en déduit l'alternance des signes et la décroissance 
des valeurs absolues. 

Nous ne considérons que les tels que 3fî(F(— w^)) < car < 0. A 
chaque U2ki on associe alors un 7^ tel que c^^ = F{—iu2k), i-^- 

7.(Ao, Al) := 1 + ^-^^^^ \F{-iU2u)\-' • 

La décroissance des valeurs absolues des F{—iu2k) montre que les 7^ sont 
ordonnés par ordre croissant. 

Nous venons de montrer que lorsque 7 varie, une paire de valeurs propres 
traverse l'axe imaginaire en chaque ^k-, et seulement en ces points-là. Lorsque 
7 tend vers 1 par valeurs supérieures, décroît vers —00, et donc les valeurs 
propres A doivent rendre F{\) de plus en plus grand en module et négatif. Or, 
le module de -F(A) est borné sur le demi-plan 3Ç(A) > (d'après l'équation 
fl3.26p ). donc pour 7 assez proche de 1, toutes les valeurs propres ont une 
partie réelle négative. C'est donc le cas pour tout 7 < 70- □ 
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En réalité, 70, 7i, • • • dépendent peu de Aq, et leurs valeurs typiques sont 
7o ^ 6,2 et 7i ^ 30. 

Lorsque 7 grandit et traverse 70 (Aq étant fixé), on s'attend à voir une 
bifurcation de Hopf (voir annexe IC.6.4|) : une orbite périodique attractive est 
créée au voisinage de l'équilibre pour 7 = 70, et attire toutes les solutions 
proches de l'équilibre (sauf l'équilibre lui-même) pour 7 > 70 proche de 70. 

4 Simulations numériques 

La mise en œuvre de simulations du modèle (12. ip a demandé un travail 
préliminaire de discrétisation et de mise au point des paramètres de simula- 
tion qui est détaillé en annexe \M Les précisions concernant le traitement des 
données sont données en annexe [B| dans l'ordre de présentation des résultats. 
Elles sont cependant indispensables pour une réelle compréhension de ceux- 
ci, car les nombreuses approximations qui ont été nécessaires ont souvent une 
réelle influence sur les résultats obtenus. 

4.1 Explorations de l'espace des paramètres 

Pour commencer, on fait varier un paramètre en gardant les autres fixes, et 
l'on observe la façon dont la dynamique stationnaire évolue. Trois paramètres 
semblent déterminants : Aq (qui introduit un effet de retard correspondant 
au temps de maturation), p (qui mesure l'importance du facteur saisonnier) 
et 7 (qui traduit l'influence de la densité sur la fécondité). Pour les autres 
paramètres, on a fixé les valeurs suivantes : 

- Ai = 2. 

- rriQ = 50 et la fécondité est prise C^. 

- eete = 0,1 (ou dans le premier cas, car le facteur saisonnier n'avait 
pas encore été régularisé). 

On pourra ainsi représenter chaque simulation par le triplet {Aq; p; 7) qui lui 
correspond. 

4.1.1 Aq = 0.18, p = 0.41, 7 variable 

Pour la première exploration, nous sommes partis des valeurs (0,18 ; 0,41 ; 
8,25) et nous avons fait varier 7. Contrairement aux simulations effectuées 
ultérieurement, le facteur saisonnier n'est pas régularisé (z.e. e = 0). Les 
résultats sont représentés sur un diagramme de bifurcation, figure [5] : pour 
chaque valeur de 7 sont représentées les valeurs de N{t) aux temps entiers 
(z.e. à la fin de l'été), en se limitant à t assez grand (on a fixé arbitrairement 
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19001 < t < 20000). Les points bleus correspondent à une même condition 
initiale (obtenue aléatoirement), que nous notons (I) (voir figure HH])- La 
partie rouge correspond à d'autres simulations, détaillées ci-après. 

Orbites périodiques attractives Pour 2 < 7 < 7,2, le comportement ob- 
servé est assez classique. On a d'abord un équilibre attractif, puis une orbite 
de période 2, et ainsi de suite avec des doublements de périodes successifs (de 
plus en plus rapprochés) au fur-et-à mesure que 7 se rapproche de la valeur 
limite 700,1 ~ 7,36 . Il s'agit de bifurcations par doublement de période (voir 
la section IC.6.2P qui ont lieu pour 70 < 71 < • • • < 7„ < • ■ ■ < 7oo,i, et qui se 
traduisent par des doublements de période successifs. 

Pour 7 > 7oo,i, le diagramme de bifurcations permet de distinguer es- 
sentiellement deux comportements. D'une part, il y a toujours des orbites 
périodiques attractives sur certains intervalles de valeurs de 7, tout comme il 
y a des «fenêtres de périodicité» dans le cas des polynômes quadratiques (voir 
section [Ô7|) . C'est en particulier le cas pour 9,00 < 7 < 9,39, 9,67 < 7 < 9,86 
et 12,23 < 7 < 16. On peut situer plus précisément ces fenêtres à l'aide d'un 
calcul de dimensions fractales. 

Dimensions fractales On peut aisément calculer une valeur approchée 
de la dimension fractale des ensembles limites correspondant aux différentes 
valeurs de 7. Le graphique obtenu est représenté figure O 

Les points où la dimension fractale est nulle correspond aux orbites pé- 
riodiques attractives, ce qui nous permet de les détecter bien plus facilement 
qu'en observant le diagramme de bifurcation. Dans le domaine intermédiaire 
(7,36 < 7 < 12,23), on constate qu'il y a alternativement des attracteurs de 
dimension environ égale à 1 et des fenêtres de périodicité. On peut raison- 
nablement penser qu'il existe un ouvert dense dans l'espace des paramètres 
pour lequel il existe une orbite périodique attractive^. C'est ouvert n'est en 
revanche certainement pas de mesure totale. 

La précision de ces calculs étant limitée par le faible nombre de points 
considérés pour chaque ensemble, il n'est pas aisé de déterminer s'il y a ou 
non réellement des attracteurs de dimension non-entière. Si oui, elle n'est pas 
très grande, certainement inférieure à 1,5, et probablement supérieure à 1, 
dans la mesure où le calcul effectué sous-estime légèrement la dimension de 
l'attracteur (en particulier, le petit nombre de points utilisés peut créer de 
nombreux «trous» correspondant à la mesure physique sur l'attracteur, et 
non à un trou réel dans sa géométrie). 

11. qui peut éventuellement coexister avec un autre attracteur. 
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= 2 : E = ; rrip = 50 ; m(N) ; 100 pas da bamps par an 
=0.18 ; p =0.+1 \ -(VRhBb\s ; 1900CEi £ 19999 




Figure 5 - Diagramme de bifurcation (0,18; 0,41; 7)2<7<i6- Pour certaines 
valeurs de 7, plusieurs attr acteur s coexistent. On les a obtenus en utilisant 
la condition (I) (en bleu) et en prolongeant aussi loin que possible (en rouge) 
certaines branches interrompues dans le diagramme bleu. Noter que les sai- 
sons sont discontinues {eete = 0). 
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Evolution de la dimension fractale :A|j-0.18;p-0.41 ;e = 



10 12 14 16 



Figure 6 - Dimension fractale des attracteurs (0,18; 0,41; 7)2<7<i6- On s'est 
limité à la condition initiale (I). 



Il y a un autre argument théorique en faveur d'une dimension fractale 
supérieure à 1 lorsqu'elle n'est pas nulle. En effet, s'il y a une orbite périodique 
hyperbolique instable, sa variété instable est contenue dans l'attracteur, qui 
doit donc avoir une dimension au moins égale à 1. Il est donc difficilement 
concevable que dans un grand domaine de l'espace des paramètres on puisse 
avoir un attracteur de dimension fractale comprise strictement entre et 1. 

Attracteurs de type Hénon En-dehors des fenêtres de périodicité appa- 
raît un comportement stationnaire non-périodique, le long d'un attracteur 
qui semble constitué de deux morceaux de courbes. Un exemple est repré- 
senté avec la figure [TJ où l'on a tracé les points {N{t), N{t + 1), N{t + 2)) 
pour chaque valeur entière de t (10002 < t < 19999). On s'intéresse alors à 
la dynamique de l'applicationEl / : {N{t), N{t + 1), N{t + 2)) ^ {N{t + 1), 
N{t + 2), N{t + 3)). Par abus de notation, on écrira au lieu de /, sans 
perdre de vue que nous ne pouvons pas visualiser directement T^. En utili- 
sant deux couleurs suivant la parité de t, on constate que chacune des deux 
parties de l'attracteur est envoyée sur l'autre. En revanche, il ne semble pas 
possible (pour cette valeur de 7) de séparer de la même façon l'attracteur en 

12. Cette application n'est pas parfaitement bien définie, le système étudié étant de 
dimension infinie, alors que la visualisation considérée est une projection de celui-ci en 
dimension 3. 
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=0.18 A, = 2 ;p = 0.4-1C E = D ; rr^ = 50 -y = 8.61 m(N) 
10001 1999B : 1.2S35£ N£â.1125 ; 100 pas datamps par an 



5.5^. ■ 

4.5^ ..■ 
4^ .. 




2^ - ■ 
1 .5 ^ . ■ 




Figure 7 - (0,18; 0,41; 8,61). Condition initiale (I). Les deux composantes 
sont invariantes par T^, qui y semble topologiquement mélangeante. Dimen- 
sion fractale estimée '■ df ^ 1,06. 



un plus grand nombre de composantes. Il semble donc que l'application 
restreinte à chacune des deux composantes de l'attracteur soit topologique- 
ment mélangeante (voir définition IC. 12p . 

Un zoom sur l'attracteur fait apparaître des structures semblables à celles 
de l'attracteur de Hénon (voir annexe IC.3.71) . Il est assez probable que la 
dynamique soit du même type, mais cela n'à pas été étudié précisément. La 
dimension fractale n'est pas clairement différente de 1, mais elle est sûrement 
sous-estimée à cause du petit nombre de points que nous avons calculé. Si elle 
s'avérait être clairement différente de 1, cela renforcerait l'hypothèse d'une 
dynamique de type Hénon. 

Il n'y a pas toujours deux composantes. On le voit sur le diagramme de bi- 
furcations notamment au voisinage des valeurs «limites» de 7 (z.e. juste après 
une zone où il y a une orbite périodique attractive). Ainsi, pour 7 = 8,62, 
on peut compter 10 composantes distinctes, et semble topologiquement 
mélangeante sur chacune d'entre elles (figure [8]). Dans les deux cas (figures 
[7] et [8]), il semble donc qu'on ait une décomposition spectrale (voir théorème 
IC.12P de l'attracteur A en un nombre fini de composantes (respectivement 2 
et 10). 

Lorsque 7 évolue!^, on observe que l'attracteur grandit petit-à-petit à 
13. voir aussi l'amination f ilm_gcLinma_100_18_200_410_00_50 1 .avi (diagramme de 
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=0.18 = Z ; p = 0.410 E = ; m^^ = 50 y = 9. 88 m(N) 
10O01 £t£ 19B9B ; 1 .233+ £ N S +.64-27 ; 100 pas detemps par an 



+.5 




Figure 8 - (0,18; 0,41; 9,89). Condition initiale (I). Chacune des compo- 
santes est invariante par T^°, qui y semble topologiquement mélangeante. La 
numérotation des composantes correspond aux valeurs de tmod. 10 . Comme 
à la figure [TJ il y a deux groupes suivant la parité de t. Dimension fractale 
estimée : df ~ 0,92. 
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l'intérieur d'un même objetll^, que l'on visualise à peu près en entier avec 
la figure [H Tant que 7 < 9,96, on remarque même que les orbites pério- 
diques attractives sont contenues dans le même objet. Lorsque l'attracteur 
est continu mais en morceaux distincts, envoyés périodiquement l'un dans 
l'autre, chaque morceau grandit et ceux-ci fusionnent petit-à petit. Lorsque 
deux composantes fusionnent, la stabilité de chacune vis-à-vis de / (com- 
posée le bon nombre de fois) semble instantanément perdue, et il y a alors 
mélange topologique à l'intérieur de chaque composante. Cette fusion des 
composantes connexes s'apparente à la cascade inverse qui suit la cascade 
harmonique directe, comme c'est le cas pour la famille quadratique réelle 
(annexe [073) ■ 

Discontinuités du diagramme Les discontinuités observées à 7 = 3,70 et 
7 = 12,23 sont plutôt surprenantes. L'hypothèse la plus plausible serait qu'il 
existe à ces valeurs de 7 deux attracteurs distincts, et la condition initiale 
(I) passe brusquement du bassin de l'un au bassin de l'autre. Pour tester 
cette hypothèse, nous avons choisi comme nouvelles conditions initiales les 
états stationnaires obtenus de part et d'autre de ces discontinuités, et nous 
avons fait varier 7 pour déterminer s'il y a effectivement coexistence de deux 
attracteurs pour certaines valeurs de 7. 

En repartant de l'équilibre obtenu avec 7 = 3,70, nous avons pu continuer 
la branche du diagramme jusqu'à 7 = 4,20. Ensuite, on retrouve l'orbite de 
période 2 déjà trouvée. Le même procédé nous permet de continuer jusqu'à 
7 = 4,23, mais nous n'avons plus pu retrouver l'équilibre ensuite en utilisant 
le même procédé. Dans l'autre sens, on prolonge le domaine où se trouve 
une orbite de période 2 attractive jusqu'à 7 = 3,60. On peut ainsi tracer un 
nouveau diagramme de bifurcations autour de ces valeurs, avec cette fois les 
deux états stationnaires (figure [9]). 

La seconde discontinuité a lieu autour de 7 = 12,23 : on observe l'attrac- 
teur pour 7 = 12,22 et une orbite attractive de période 4 pour 7 = 12,23. En 
repartant de l'état final obtenu pour 7 = 12,22, on observe pour les valeurs 
supérieures de 7 (au moins jusqu'à 16) un comportement similaire à ce qu'on 
constatait pour 7 < 12,22, c'est-à-dire le même attracteur, avec parfois des 
orbites périodiques attractives (mais dans un domaine de valeurs de dif- 
férent de l'orbite de période 4). Inversement, l'orbite 4-périodique attractive 
persiste jusqu'à 7 ~ 10,078. 

Les attracteurs ainsi détectés sont représentés en rouge sur la figure [5l 
Dans ces deux cas, plusieurs attracteurs coexistent, mais il y a toujours des 

bifurcation en quatre dimensions). 
14. A de légères déformations près, notamment un changement de taille. 
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A^=2;e^0;m^ = 50;m(N)C\l00pas de temps paran 
= 0.18 :p - 0.41 ; Y variable ; 19001<t< 20000 



Figure 9 - Diagramme de bifurcation (0,18; 0,41; 7)3,5o<7<4,30- Deux attrac- 
teurs coexistent pour certaines valeurs de 7 : un point fixe (au centre) et une 
orbite de période 2 (en haut et en bas). Avec la condition (I), on passe de 
l'un à l'autre pour 3,70 < 7 < 3,75. 



discontinuités dans le diagramme. Pour certaines valeurs de 7 (envion 3,60 
et 10,078), des orbites périodiques deviennent attractives. A l'inverse, pour 
7 4,23, l'équilibre devient instable. Il pourrait donc s'agir d'une bifurcation 
du type de /(+i),r (voir annexe [C6.2[ figure [30]) • H est également possible que 
le diagramme soit en réalité continu, mais que les bassins d'attraction des 
orbites périodiques attractives soient trop réduits pour que l'on puisse les 
atteindre par des simulations, avec la méthode que nous avons employée ici. 
Un petit travail théorique serait nécessaire pour éclairer ce point. 

4.1.2 Aq = 0,18, p = 0,30, 7 variable 

Suite à une rapide exploration en faisant varier p, et au vu de l'intérêt des 
valeurs (0,18; 0,30; 8,25) (voir figure fT2|) . nous avons effectué une deuxième 
exploration à 7 variable, autour de ces nouvelles valeurs. Le diagramme de 
bifurcation ainsi obtenu est représenté figure [101 H est très semblable au 
diagramme 13 mise à part l'absence de deux composantes bien distinctes pour 
la plupart des valeurs de 7, et le faible nombre de fenêtres de périodicité. 
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= 2 ; E = 0.1 ; = 50 : m(N) : 100 pas datemiis par an 
Ajj= O.IB :p = 0.3 Invariable ; 1 900 l£t 5.20000 





Figure 10 - Diagramme de bifurcation (0,18; 0,30; 7)2<7<i5- Les saisons sont 
continues {eete = 0,1). Condition initiale (I). 
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= 2 ;e = 0.1 ; = 50 ; m[N) ; 100 pas da temps par an 
fl^ = Û.ia ; p variable iT= 6.25 ; 190Q1S t £ 20000 

T 1 1 1 1 1 r 




0.06 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 



P 

Figure 11 - Diagramme de bifurcation (0,18; p; 8,25)o<p<o,5- Les saisons sont 
continues {eete = 0,1). Condition initiale (II). 



4.1.3 Ao = 0,18, p variable, 7 = 8,25 

Le facteur saisonnier semble jouer un rôle déterminant dans la dynamique 
du système (le modèle non-saisonnier est particulièrement simple, alors que 
pour des valeurs de p plus proches de la réalité, on observe des comportements 
bien plus complexes, e.g. figure [T2|) . Le diagramme de bifurcation obtenu en 
faisant varier p est représenté figure [TTl Une autre condition initiale, notée 
(II), a été utilisée pour ces simulations (voir figure lïïU]) . On observe comme 
précédemment une discontinuité dans le diagramme, pour p proche de 0,1, 
mais nous n'avons pas essayé de prolonger les deux branches interrompues. 

Bifurcation de Hopf On constate sur ce diagramme un phénomène que 
nous n'avions pas trouvé sur les diagrammes précédents : une bifurcation de 
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A^j= 0.1B00 A, =2 :p= 0.300 e =0.1 ; m^=5(l ï= B.S mfN) 
1 0002 Sts 19998 ; 1.1933 S N£ 6.0647 ; 1O0 pas de tampB par an ; ^ non-enlier 




Figure 12 - Un comportement complexe : (0,18; 0,30; 8,25). Les saisons sont 
continues {este = 0,1). Condition initiale (II). Il semble que soit mélan- 
geante sur cet attracteur. Dimension fractale estimée : df ~ 1,19. 
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A^j= O.ISA, = 2 ;p = 0.160 e =0.1 ; m^^ = 50 ï = B.25 m(N) C"' 
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Figure 13 - Deux «cycles» : (0,18; 0,16; 8,25). Condition initiale (II). Chaque 
composante est invariante par T^. Dimension fractale estimée : dj ^ 0,99. 
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Hopf pour p ^ 0,152. On passe en effet d'une orbite attractive de période 2 à 
deux «cycles» attractifs stables, qui persistent jusquà p ~ 0,165 (la figure [T3] 
en représente un exemple). 

Cet attracteur est de dimension fractale 1 et chaque lacet est parfaitement 
connexe. Il n'est en revanche pas totalement certain que soit bien topolo- 
giquement mélangeante sur chacun, bien que l'on n'ait vu aucune périodicité 
évidente. La dynamique de / sur ces cycles n'est pas forcément simple {i.e. 
topologiquement conjuguée à une rotation), et le cycle ne coïncide peut-être 
pas exactement avec l'attracteur. Répondre à ces questions demanderait une 
étude plus poussée. 

4.1.4 Aq variable, p = 0,3, 7 = 8,25 

Pour cette dernière exploration préliminaire, une difficulté supplémentaire 
a été de trouver une façon de faire varier la valeur de Aq plus finement que le 
pas de discrétisation. Celle-ci a conduit à utiliser la méthode évoquée à la fin 
du paragraphe lA.ll Les résultats obtenus sont représentés dans le diagramme 
de bifurcations de la figure [Hl 

Décomposition spectrale L'utilisation de 2 couleurs de visualisation montre 
que la 2-périodicité de la figure [TSl ne se retrouve pas à la figure fT2l II y ainsi 
initialement deux composantes connexes bien distinctes, l'une étant l'image 
de l'autre par l'application T^. Sur chaque composante, semble mélan- 
geante. On a ainsi une décomposition spectrale avec 2 composantes (voir 
théorème IC.12|) . Lorsque celles-ci fusionnent, on perd cette 2-périodicité et 
devient topologiquement mélangeante. On retrouve le même comporte- 
ment que précédemment avec le diagramme (0,18; 0,41; 7). 

4.2 Étude du cas (0,15; 0,30; 8,25) 

Essayons de comprendre la dynamique de la figure [151 Nous avons vu 
qu'il y a deux composantes connexes distinctes, il suffit donc de considérer 
l'une des deux pour comprendre la dynamique de T^. Elle est représentée 
figure [ini 

4.2.1 Visualisation en dimension 3 

La représentation que nous avons choisie (arbitrairement) est-elle cor- 
recte ? Cette question est fondamentale : nous projetons en effet un objet de 
dimension infinie!^ dans un espace de dimension 3. 

15. après discrétisation, on se limite à une dimension finie très grande, ici 201. 
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= £ ; E = 0.1 ; = 50 ; m(N) : 100 pas da temps par an 
variable ; p = D.3 ; y = 8.25 ; 19001at £ 20000 




t 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 



Figure 14 - Diagramme de bifurcation [Aq] 0,30; 8,25)o<^o<o,4- Condition 
initiale (11). 
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5.5^ ■• 




Figure 15 - Un comportement complexe : (0,15; 0,30; 8,25). Condition ini- 
tiale (II). Les deux composantes sont invariantes par T^, qui semble topolo- 
giquement mélangeante sur chacune. 
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Figure 16 - (0,15; 0,30; 8,25) Condition initiale (II). On n'a conservé qu'une 
composante connexe. 
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Test de l'injectivité de la projection L (R 



» L'"'(R'') 




1000 ~ 



30 40 50 

j = n° sur la carte 



Figure 17 - Injectivité de la projection : norme ôt = 0. 



Injectivité de la projection Pour tenter d'y répondre, nous pouvons 
évaluer la qualité de la «projection» vr : M^^-*^ — > R^, 0:201 (^) = {N{t + 
fc/100))fc=o...2oo ^ {N{t), N{t + 1), N{t + 2))= xsit) où N{t) désigne la popu- 
lation mature à l'instant t, en différents points de l'attracteur. Nous voulons 
nous assurer que des points proches dans sont également proches dans 
R^°\ c'est-à-dire majorer sup^-^^/^pj ^ ^ ^ 1 "xi (t) -x" (f T 3^°^ pour différents choix de 
normes (L^, ou 

Le résultat, représenté à la figure [T71 montre que cette quantité est rai- 
sonnablement bornée. Avec les normes ou L^, le résultat semble un petit 
peu meilleur, mais reste du même ordre de grandeur. Une zone de l'attracteur 
semble en revanche être un peu moins bien représentée par cette projection, il 
s'agit du point 24 (et plus généralement des points 21 à 30). En se reportant 
à la figure |25] où sont localisés ces points (voir section 14.2.31) . on constate 
qu'il s'agit de la zone de pli. Une vue rapprochée sur cette zone de l'attrac- 
teur montre en effet des filaments entrelacés, et certains rapprochements de 
filaments semblent dûs à la projection. 

Choix de l'origine des temps Nous l'avons arbitrairement fixée à la 
fin de l'été, mais ce choix est-il judicieux ? Nous avons donc fait les mêmes 
calculs que précédemment en décalant l'origine des temps. Il semble que 
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Injectivité : choix de l'origine des temps 

60 1 , , , , , , , , , ^6 




Figure 18 - Injectivité de la projection : norme L°°, St variable. 



l'instant choisi initialement n'est pas mauvais. La «meilleure» origine semble 
se situer autour de ôt = 0,4, mais la différence avec St = n'est pas ffagrante 
(figure [HD. 

On peut expliquer ces résultats en observant l'évolution en temps continu 
de N{t). En effet, la population mature atteint tous les deux ans — un 
peu après le milieu de l'été — un maximum élevé, suivi d'une chute brutale 
d'effectif. La valeur ôt = —0,4 correspond à l'instant du pic de population, 
qui est suivi d'une simple diminution linéaire de N{t) (dûe à la mortalité 
naturelle, en l'absence de naissances), si bien que les instants qui suivent 
sont encore des origines des temps de bonne qualité. 

Il est également intéressant, en vue de comprendre la dynamique en temps 
continu du système, de visualiser l'évolution de l'attracteur tridimensionnelle 
lorsque l'on fait varier l'origine des temps ôt dans l'intervalle [—1; 1], la valeur 
correspondant à la fin de l'été. La figure [19] en donne un exemple, pour 
ôt = -0,4. 

Echelle logarithmique Une autre piste possible est de visualiser la pro- 
jection tridimensionnelle de l'attracteur suivant une échelle logarithmique, 
c'est-à-dire de considérer les points (log A^(t), log A^(t + 1), log A^(t + 2)) pour 

16. l'animation filin_delta.avi permet mie bonne compréhension de la façon dont 
l'attracteur se déforme, pom' passer d'une composante à l'autre quand ôt varie de —1 à 
ou de à 1. 
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Figure 19 - Visualisation de l'attracteur (0,15; 0,30; 8,25) avec une origine 
des temps 5 = —0,4. 
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Figure 20 - Visualisation de l'attracteur (0,15; 0,30; 8,25) avec une échelle 
logarithmique. 



t entier grand. 

Au vu de la figure [201 l'apport d'un tel changement d'échelle n'est pas 
évident. L'attracteur est légèrement déformé, mais garde le même aspect, et 
certaines zones semblent toujours aussi «emmêlées». Une évaluation quantita- 
tive de la qualité de cette nouvelle projection, comme effectué précédemment, 
confirme l'aspect visuel : il n'y a pas de gain significatif. 

4.2.2 Géométrie de l'attracteur 

La figure 122] représente de façon simplifiée la géométrie de l'attracteur de 
la figure [161 ^n distinguant neuf régions principales. Celles-ci sont nommées 
en fonction de leur forme et de leur position dans l'attracteur, identifié au 
corps d'un animal dont la tête serait située à droite (chevelure-cou-pli-pointe) 
et la queue à gauche. 
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Injectivité : choix de l'origine des temps 

400 I \ \ \ 1 \ \ \ \ \ ^8 




q\ \ \ \ 1 \ \ \ \ \ ^4 

-0.5 -0.4 -0.3 -0.2 -0.1 0.1 0.2 0.3 0.4 0.5 



Figure 21 - Qualité de la projection en échelle logarithmique. 



Quatre phénomènes principaux semblent pouvoir expliquer la complexité 
de l'attracteur que nous observons. Les deux premiers étaient déjà présents 
dans le solénoïde : un très fort pincement et un étirement. Le troisième est 
également présent dans l'attracteur de Hénon : un pli (il y en a peut-être 
plusieurs ici). Le quatrième semble nouveau, et ressemble à un ou plusieurs 
«embranchements». Avec la numérotation introduite dans la section précé- 
dente, on peut en situer trois : 67-80, 63-55, 58-74. Cela ne signifie pas pour 
autant que ces embranchements sont distincts. 

Localement, l'attracteur ressemble au produit d'une droite et d'un en- 
semble de Cantor (figure [23|) . sauf en certains points où l'on observe des 
«pointes» (figure [23]) . Par ces aspects, il ressemble beaucoup à l'attracteur 
de HénonPU 

4.2.3 Dynamique sur l'attracteur 

Considérons l'application T^. Comment agit-elle sur les points de l'at- 
tracteur ? Pour essayer de le comprendre, 80 points ont été choisis!^ sur 
l'attracteur, numérotés de 1 à 80, comme représenté sur la figure [25] (des 
vues plus rapprochées sont en annexe ]D]) . La position des images directes et 
réciproques de ces 80 points est indiquée dans le tableau [T] 

17. voir annexe [Ç. 3. 71 

18. Ce choix a été fait arbitrairement, en essayant de répartir ces points uniformément 
suivant la mesure de Hausdorff sur l'attracteur, et non la mesure physique. 
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Madàl9 géomèlriq U9 simplifié' d'attractaur 



5.5-, 




Figure 22 - Représentation grossière de la géométrie de l'attracteur 
(0,15; 0,30; 8,25). 



Figure 23 - Zoom sur un filament, au voisinage du point 23. 
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26-27 (9^+) ^ 13 ^ 69 




55- 


-62 ^ 
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^ 38 




27-15 ^ 14 46- 


-34 




53 ^ 


5 


^ 39- 


-40 


39 ^ 15 ^ 28- 


-16 




62 ^ 


6 


^ 39- 


-40 


32-23 ^ 16 ^ 28- 


-27 


62- 


-64 ^ 


7 


^ 39- 


-40 


21 ^ 17 ^ 35 






52 ^ 


8 


^ 30 




24 ^ 18 ^ 69- 


-69 


51- 


-75 ^ 


9 


^ 40 




29 ^ 19 ^ 46- 


-34 


59- 


-60 10 ^ 68 




40 ^ 20 ^ 26- 


-27 





46 ^ 21 ^ 17 




45- 


-46 ^ 31 ^ 24 




43- 


-43 ^ 41 ^ 53- 


-55 




34 ^ 22 ^ 18- 


-13 




37 32 ^ 16- 


-28 




45 ^ 42 ^ 62- 


-55 




50 ^ 23 ^ 18- 


-17 


46- 


-45 ^ 33 ^ 15- 


-17 




49 ^ 43 ^ 41- 


-42 (13+) 




31 ^ 24 ^ 18 




12- 


-13 34 22 




13- 


-14 ^ 44 ^ 62- 


-64 


39 


-38 ^ 25 ^ 20- 


-19 




17 ^ 35 ^ 11 
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-18 ^ 45 ^ 42 






12 ^ 26 ^ 13- 


-12 


14- 


-15 ^ 36 ^ 12 
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20- 
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-76 ^ 47 ^ 61- 


-63 


15- 


-12 28 ^ 14- 
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04 ^ 38 ^ 24- 
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-06 ^ 48 ^ 63- 


-41 
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-06 29 19 
08 30 ^ 18- 
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03 39 15 
09 ^ 40 ^ 20 






01 49 43 

02 ^ 50 ^ 23 





41- 


-63 ^ 51 ^ 09- 


-08 




68 ^ 61 ^ 03- 


-07 


63- 


-64 ^ 71 ^ 80 




41- 


-63 ^ 52 ^ 08 




70- 


-45 ^ 62 ^ 06 






65 ^ 72 ^ 79 




41- 


-63 ^ 53 ^ 05 




43- 


-41 ^ 63 ^ 07- 


-03 


63- 


-41 ^ 73 ^ 78 






69 -> 54 ^ 01- 


-05 


44- 


-46 ^ 64 ^ 75 






66 ^ 74 ^ 77 




41-42 (1 


3+) ^ 55 ^ 04- 


-02 


77- 


-78 ^ 65 ^ 72 






64 ^ 75 ^ 76 




67- 


-80 ^ 56 ^ 08- 


-03 


77- 


-76 ^ 66 -> 74 






75 ^ 76 ^ 47- 


-67 




80 ^ 57 ^ 09- 


-08 


71- 


-07 ^ 67 ^ 56- 


-61 




74 ^ 77 ^ 65- 


-66 


66- 


-67 ^ 58 09- 


-76 




10 68 ^ 61 






73 ^ 78 ^ 60 






79 ^ 59 ^ 10- 


-09 




13 69 -> 54 






72 ^ 79 ^ 59 






78 ^ 60 ^ 10- 


-72 


18-17 (23+) ^ 70 ^ 62- 


-54 




71 ^ 80 ^ 57 





Table 1 - Dynamique sous des 80 points sur l'attracteur. On a noté 1^+ 
pour T®(1) ; 54-55 pour un point situé entre 54 et 55, plus proche de 54; 
76-76 pour un point proche de 76 mais différent de 76. 
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Figure 24 - Zooms successifs sur une pointe, au voisinage du point 54. 



On peut penser qu'il y a un point fixe dans la région du point 15, une 
orbite de période 2 dans la région 42-31 et une orbite de période 3 dans le 
«triangle» 66-74-77. Nous étudierons le point fixe plus en détails dans une 
des sections suivantes. Partant de ces trois orbites remarquables, on peut 
essayer de comprendre schématiquement la dynamique interne à l'attracteur. 

Sortie de la «tête» L'équilibre (entre pli et pointe) est répulsif avec ex- 
pansion négative (la valeur propre dominante est négative), tandis que l'or- 
bite de période 2 (dans le cou, entre les points 42 et 31) est répulsive positive. 
Dans cette région de la tête et de la pointe arrière, on peut schématiser la 
dynamique de la façon suivante : 

- Pointe arrière — creux-cou-pli-pointe (sans inversion haut-bas). 

- Chevelure — )■ cou-pli-pointe (sans inversion). 

- Bas du cou (sous l'orbite de période 2 = 42-31) — j- plat-creux-bas du 
cou (sans inversion). 

~ Haut du cou (au-dessus de la période 2) — )■ haut du cou-pli (sans in- 
version) . 

- Pli — )■ bas de la pointe (avec inversion). 

- Haut de la pointe — )■ bas du pli (avec inversion). 

Ainsi, si l'on part d'assez haut, on arrive en bas de la pointe puis en bas 
de la chevelure (69) et enfin dans le creux (54). Si l'on part plus bas, on arrive 
directement dans le creux, sans être passé par la pointe. 

Sortie du «corps» Dans le «corps» (queue, triangle, plat, creux), l'orbite 
de période 3 (aux milieux des sommets du triangle) joue un rôle majeur. Elle 
est répulsive positive, et agit sur le triangle comme une rotation d'angle 27r/3 



42 




.71: 



P5 



<Ç4 



Figure 25 - Position des 80 points choisis sur l'attracteur. 
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(75-72 —7- 76-79 — )■ 67-59) combinée avec un peu d'expansion. On sort ainsi 
du triangle soit par le bas de la pointe arrière, soit par le plat ou le creux, 
soit par la queue. La dynamique se schématise alors ainsi : 

- Bas de la pointe arrière — ?• plat-creux. 

- Creux — )• droite de la queue (l-7)-côté du triangle (75-73) (avec inver- 
sion droite-gauche). 

- Plat — > gauche de la queue (sans inversion haut-bas). 

- Queue — )■ pointe arrière. 

On sort ainsi du corps, pour y revenir rapidement (si l'on arrive trop bas 
dans la pointe arrière, soit près du triangle, soit à la pointe 49), ou (le plus 
souvent) après un passage dans le pli et éventuellement la pointe. 

Ce bref aperçu de la dynamique nous permet de comprendre comment 
s'instaure le mélange et le chaos de la dynamique sur l'attracteur. 

4.2.4 Visualisation en temps continu 

Il est intéressant, du point de vue mathématique comme du point de vue 
biologique, de mettre en relation les différentes régions de l'attracteur (en 
dimension 3) avec la dynamique en temps continu dans ces régions. La forme 
de N{t) pour t G [to — 5; to + 5] peut ainsi être mise en correspondance avec 
les points de l'attracteur et leurs régions de provenance et de destination par 
T^. On a représenté ces données pour les 80 points de la «carte» en annexe iDl 

Aspect général La dynamique en temps continu garde en général un as- 
pect identique dans toutes les zones de l'attracteur : un pic chaque année 
pour t ~ 0,6 mod. 1 (croissance très rapide suivie d'une décroissance linéaire 
moins brusque), élevé les années impaires (entre 2,5 et 6,5), plus faibles voire 
inexistant les années paires (entre 1 et 2). Juste avant les maxima se trouvent 
des minima locaux plus ou moins bas : ceux-ci sont toujours légèrement su- 
périeurs à 1 (mais inférieurs à 1,5) avant un maximum faible, mais pouvant 
atteindre 0,5 avant un maximum élevé. Il y a ainsi une pseudo-périodicité de 
2 ans, avec une très forte variation d'amplitude (les maxima étant à 6,5, les 
minima entre 0,5 et 1). On retrouve en partie l'aspect général de la figure H] 
(alternance de maxima et de minima, avec un facteur allant jusqu'à 8 entre 
les deux). 

Ce comportement s'explique par l'alternance entre une explosion de la 
population dûe à la très forte fécondité, qui est ainsi suivie d'une chute li- 
néaire de la population mature (dûe à la mortalité naturelle, en l'absence 
quasi-totale de naissances). 
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Différences entre régions Les régions de l'attracteur se différencient par 
l'amplitude des pics (il y a un facteur 2 entre les amplitudes possibles des 
maxima élevés) et des creux (inférieurs ou supérieurs à 1), ainsi que par l'am- 
plitude relative des maxima secondaires (inexistants ou bien valant jusqu'à 
2). 

Les pics faibles se trouvent au niveau de la queue (lorsqu'il se produit à 
t = 0) ou du plat (lorsqu'il se produit k t = 2). 

Les pics très élevés (A^ > 6) se trouvent dans le pli, la pointe, le haut de la 
pointe arrière et le cou, avec des maxima secondaires quasi- inexistant s et des 
minima inférieurs à L Le deuxième pic est légèrement inférieur au premier 
dans la pointe, mais le rapport s'équilibre quand on se rapproche du pli. 

On peut encore affiner cette analyse en s'aidant des figures placées en 
annexe [Dl qui permettent de faire la différence entre fflaments dans une 
région à l'aide de l'évolution en temps continu dans le passé ou dans le futur 
proche. 



4.2.5 Dimension fractale de l'attracteur 

Pour l'évaluer, nous avons calculé pour différentes valeurs de r le nombre 
de boîtes de côté r (et dont les points ont des coordonnées qui sont des 
multiples entiers de r) contenant des points de l'attracteur (avec 100 000 
points). En notant N{r) ce nombre de points, on a alors tracé logiQ{N{r)) 
en fonction de log]^Q(r). Tant que N{r) est assez petit devant 100 000 et 
assez grand devant 1, les points obtenus sont presque alignés, et la pente 
(en valeur absolue) de la droite de régression est une bonne estimation de la 
dimension fractale de l'attracteur. Le résultat est représenté à la figure [221 
On obtient donc une dimension de l'ordre de 1,33. Cette valeur correspond 
bien à l'impression visuelle que l'on a : localement, l'attracteur semble être 
le produit d'une droite et d'un ensemble de Cantor de dimension proche de 
1/3, soit une dimension fractale d'environ 4/3 (voir figure [23]) . Il faut bien sûr 
prendre ce résultats avec beaucoup de précautions, dans la mesure où cette 
régression est faite dans la zone qui semble — visuellement — présenter une 



«bonne pente» (voir annexe \B.6\i . 

Une des conséquences de cette évaluation est la justification a posteriori 
de la possibilité d'utiliser trois dimensions seulement pour visualiser l'at- 
tracteur. En effet, le théorème de Whitney IC.20 ^ affirme qu'un objet de 



dimension d peut être visualisé avec dimensions pourvu que N > 2d. L'at- 
tracteur ayant une dimension strictement comprise entre 1 et 1,5, l'utilisation 
de 3 dimensions pour le visualiser semble raisonnable!^. 



19. voir annexe [C. 101 

20. mais bien sûr, comme nous avons calculé la dimension fractale de la projection 
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4.2.6 Sensibilité aux conditions initiales 

Pour l'évaluer, il est intéressant de regarder la dynamique future et passée 
d'une petite boule centrée sur un point de l'attracteur. 

Les figures [271 et 128) permettent de juger du résultat au voisinage de l'équi- 
libre. En se référant à l'exemple du solénoïdeEïl, on peut tenter d'interpréter 
la forme des courbes obtenues. 

Tout d'abord, il y a clairement une très forte sensibilité aux conditions 
initiales, en tout point de l'attracteur, aussi bien dans le passé que dans le 
futur. Dans le cas de la figure |271 on constate ainsi des écarts de l'ordre de 4 
en moins de 15 ans, aussi bien dans le passé que dans le futur, alors que les 
courbes étaient initialement séparées de moins de 0,04. 

Les dynamiques futures des différents points se répartissent de façon à 
peu près homogène, au moins au cours des 10 premières années. Ceux-ci 
se séparent en effet selon leur répartition initiale dans la direction instable, 
puisque la direction stable est contractée dans le futur. L'homogénéité dans le 
futur traduit une répartition à peu près uniforme des points d'une orbite dans 
la direction instable. Cette propriété est à rapprocher du cas du solénoïde, 
où la mesure physique (qui donne la répartition des points d'une orbite sur 
l'attracteur) possède une densité par rapport à la mesure de Lebesgue dans 
la direction instable. 

La dynamique passée semble bien différente, les différentes courbes se 
séparant «par paquets», et non plus de façon homogène. Ainsi, en à peine 3 
ans, on observe déjà une différence de 1 pour l'une des courbes, tandis que 
la plupart des points ont une orbite toujours très proche de l'équilibre. Ceci 
nous donne des informations sur la mesure physique dans la direction stable, 
puisque la direction instable est contractée dans le passé. Ainsi, comme dans 
le cas du solénoïde, il semble que la mesure physique possède une densité par 
rapport à la mesure de Hausdorff sur un ensemble de Cantor de dimension 
fractale 0,3. 

Les résultats observés ici sont cependant moins clairs que dans le cas du 
solénoïde, pour |t| > 10. Ceci est sans doute dû au passage des orbites dans 
un pli, phénomène qui ne se produit pas dans le cas du solénoïde. 

On observe le même type de résultat en de nombreux autres points de 
l'attracteur, d'autant plus nettement que l'on n'est pas au voisinage d'une 
«pointe». 

tridimensionnelle de l'attracteur, nous n'avons pas la dimension fractale de l'attracteur 
lui-même mais une légère sous-estimation de celle-ci. 
21. voir annexe [Ç. 3. 61 
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Figure 27 - Dynamique passée et future de points au voisinage de l'équilibre. 
Les points ont été choisis pour leur proximité de l'équilibre dans l'intervalle 
de temps [0; 2]. 
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Figure 28 - Position des points au voisinage de l'équilibre : l'ensemble des 
points de la boule considérée est représenté en vert. 
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Figure 29 - Position de l'équilibre, et de sa «préimage». 



4.2.7 Point fixe, variété instable 

Un point fixe (instable, bien sûr) a été repéré sur l'attracteur. La figure [29] 
donne sa localisation approximative en dimension 3, et la fonction continue 
correspondante est représentée figure [5U1 

Différentielle de à l'équilibre Celle-ci nous fournit beaucoup de ren- 
seignements sur le système dynamique. On peut calculer aisément ses valeurs 
propres. Une seule est de module strictement supérieur à 1 et correspond à 
l'expansion dans la direction instable : Ai ~ —2,29. Le vecteur propre associé 
est tracé figure [SH La seconde plus grande valeur propre (la précision de 
ce calcul est faible) est A2 ~ 0,043, et le vecteur propre associé est repré- 
senté figure [521 Les modules des valeurs propres suivantes décroissent ensuite 
rapidement, comme le montre la figure [331 

Le point fixe est donc hyperbolique, puisqu'il n'a aucune valeur propre 
de module proche de L De plus, on constate que l'expansion est assez rai- 
sonnable (elle est sans doute un peu plus forte dans certaines zones, un peu 
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L'équilibre : (N(t),N{t+1),N(t+2))-(4.5381, 1.2352,4. 5381) 



6 - 




0.5 1 1.5 2 2.5 3 3.5 4 

t 

Figure 30 - L'équilibre : temps continu. 



Vecteur propre de la différentielle à l'équilibre ; ^ -2.2942 ; e = 1 e-006 



z: 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 



Figure 31 - Différentielle à l'équilibre : premier vecteur propre. 
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Vecteur propre de la différentielle à l'équilibre = 0.043251 ; e = 1e-006 

P| . 1 1 1 1 1 1 



1.5 - 



0.5 - 



- 



-0.5 - 



-1.5 - 




-2' ' ^ ^ ^ ^ ^ ^ ^ ^ ' 

0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 



Figure 32 - Différentielle à l'équilibre : second vecteur propre. 



Différentielle à i'équiiibre : décroissance des modules des valeurs propres 



100 150 200 

j 



Figure 33 - Différentielle à l'équilibre : décroissance des valeurs propres. 
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Variété instable : f ^(W'^ix )) 



5- .■ 




Figure 34 - Variété instable. 



moins dans d'autres, mais reste de cet ordre de grandeur), tandis que la 
contraction est beaucoup plus forte. Cela nous donne un argument supplé- 
mentaire pour penser que 3 dimensions suffisent à représenter l'attracteur : 
les valeurs propres suivantes ayant un module encore plus petit, l'attrac- 
teur n'est vraiment étendu que dans 2 ou 3 dimensions, les autres étant peu 
importantes. 

Variété instable On peut déterminer la variété instable en regardant les 
images par / = d'un segment situé dans la direction instable au voisinage 
de l'origine (figure . 

Il est intéressant, d'un point de vue dynamique, de visualiser comment 
celle-ci se déploie à l'intérieur de l'attracteur, à la fois d'un point de vue 
dynamique et d'un point de vue géométrique. En effet, la variété instable est 
une courbe continue, ce qui nous donne une idée plus précise de la géométrie 
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N(2t+2) 



N(2t+1) 



N(2t+2) 



N(2t+1) 



fa) n = 12 



Variété instable : f^^lWtx^^)) 



(b) n = 15 

Variété instable :f"(W''{x 




N(2t+2) 



N(2t+1) 



N(21+2) N(2t+1) 

(d) n = 17 



(c) n = 16 

Figure 35 - Déploiement de la variété instable à l'équilibre f^{W'^{Xeq)). 



de l'attracteur que lorsque nous ne disposons que d'un nuage de points. Ceci 
est fait avec une animation^ dont la figure |25] donne des extraits. 

On peut décrire f"'{W^{Xeq)) pour les valeurs successives de n de la façon 
suivante : 

- 1 < n < 11 : une petite portion rectiligne autour de l'équilibre. 

- n = 12 : un filament plié, dans la direction de la pointe intermédiaire 
située dans le cou (43) (figure [35k) . 

- n = 13 : le filament est un peu étendu et atteint le bas de la pointe. 

- n = 14 : le bas de la chevelure est atteint, ainsi que le cou. 

- n = 15 : le creux (avec la pointe 54) est atteint. Les filaments sont 
dédoublés (figure [35b). 

- = 16 : un côté du triangle (71-75) et la droite de la queue sont 

22. consulter les fichiers var_u_0 . 1_2_1_18 . avi et var_u_0 . 1_2_18 . avi. 
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atteints. Les filaments se dédoublent ailleurs (figure 15^) . 

- n = 17 : un autre côté du triangle (76-80) et la pointe arrière sont 
atteints (figure [35li) . 

- n = 18 : le dernier côté du triangle (59-66) est atteint. C'est la dernière 
région importante à être touchée. Notons tout de même que le filament 
58-74 (c'est l'un des embranchements) n'est pas encore atteint. 

4.2.8 Formation du pli 

Une des caractéristiques du système est l'existence d'un pli[^. Sous l'ac- 
tion de T, la région quasi-rectiligne de la figure [5^ se courbe progressivement 
pour arriver dans la région «pliée» de la figure \^Eb . 

Approche géométrique On peut visualiser la formation de ce pli en cal- 
culant la courbure au niveau du pli (sur un même filament de l'attracteur) 
à différents instants. Le résultats est reproduit figure [271 Hormis quelques 
irrégularités, il se forme clairement un pli pour 1,5 < t < 1,6, et celui-ci 
s'accentue fortement pour former un pli très marqué à t = 4. L'évolution du 
maximum de courbure est reportée figure [38k. Cette étude permet également 
de localiser très précisément l'endroit plié, en notant à chaque instant la posi- 
tion du maximum de courbure sur le segment. On constate avec la figure [55b 
que celui-ci est situé au point 690 du segment initialement choisi. 

Étude des discontinuités de la courbure Les différentes figures mon- 
trant la formation du pli présentent cependant quelques irrégularités. Ainsi, 
on peut se demander à quoi sont dûes les discontinuités de la courbure fi- 
gure [371 Au vu de l'imprécision de la méthode de localisation du pli, il est 
quasiment certain que la ligne brisée considérée saute d'un filament à l'autre. 
Il est également possible que le manque de régularité des fonctions du mo- 
dèle (qui sont et non C^) engendre des ruptures de courbure au niveau 
de l'attracteur lui-même. 

Pour tester cette seconde hypothèse, nous avons tenté d'utiliser la variété 
instable globale de l'équilibre, qui devrait nous assurer que l'on considère 
un seul filament. Les résultats obtenus sont représentés en annexe [Dl avec 
les figures [H21 à [SU H y a toujours des discontinuités dans la courbure le 
long du filament. On peut sans aucun doute attribuer le pic de courbure 
aux environs de j = 500 à un saut d'un filament à un autre (figure [83]) . Il 

23. Il y en a peut-être plusieurs, mais il n'est pas évident de distinguer si deux régions 
pliées sont ou non indépendantes. Cela demanderait une étude plus approfondie. 

24. Cette figure est extraite d'une animation donnant plus d'informations sur la forma- 
tion géométrique du pli. 
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(b) r2(xi),...,T2(xiooo) 

Figure 36 - Localisation du pli et de sa préimage. 
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Courbure du pli (t=0) : Maximum pour j=9S 



Courbure du pli (1=1 .5) ; Maximum pour j=49 
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j 

(a) t = 

CourburG du pli (1=2) ; Maximum pour j=680 




1000 100 200 300 400 500 BOO 700 



(b) t = 1,5 

CourburG du pli (t=2.5) ; Maximum pour j=1 25 
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(c) t = 2 



Courbure du pli (1=3.3) ; Maximum pour j=677 




(d) t = 2,5 



Courbure du pli (1=4) ; Maximum pour j=683 
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i j 

(e) t = 3,3 (f) t = 4 

Figure 37 - Formation du pli : courbure. 
57 



Formation du pli sur 4 ans 
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(a) Valeur maximale 

Localisation du pii au cours de sa formation 
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(b) Localisation du maximum 
Figure 38 - Formation du pli : maximum de courbure. 
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semble très difficile d'extraire avec suffisamment de précision un seul filament 
de l'attracteur dans la zone de pli. Cette tâche est peut-être simplement 
rendue impossible par l'aspect fractal de celui-ci. En revanche, la similitude 
des figures [371 et [831 conffime nos observations sur le processus de pliage. La 
figure [8^ est d'ailleurs plus simple à interpréter, puisque l'on a uniquement 
le pic artificiel à j ~ 500 et le pic réel à j ~ 875. 

La figure [SSlprésente également des irrégularités facilement interprétables. 
En effet, la courbure ne croît pas toujours avec le temps, et le maximum de 
courbure se déplace par moments, même après la formation du pli à t = 1,6. 
Ces irrégularités peuvent s'expliquer par le non-uniformité du pincement et de 
l'étirement dans la dynamique. Ainsi, le pincement fait augmenter fortement 
la courbure, et l'étirement entraîne au contraire sa diminution. Sur deux ans, 
cela se traduit par une forte augmentation de la courbure, mais en temps 
continu, il est des périodes où la tendance s'inverse. Un argument en faveur 
de cette hypothèse est que ces périodes coïncident avec les instants des pics 
de population, c'est-à-dire les instants d'expansion maximale (voir figures [39] 
et [12]). 

Approche dynamique On peut aborder la question du pli d'un point de 
vue dynamique et non plus purement géométrique en observant l'évolution de 
la population en temps continu (figure [5^ . Les points Xj du filament non-plié 
dans correspondent à l'intervalle [—2; 0]. 

Il est également intéressant de considérer la différentielle de le long du 
pli pour mieux comprendre ce qui se passe dans cette région de l'attracteur. 
Les vecteurs tangents (figure nous indiqueront alors comment plie le 
filament lui-même. 

Filament non encore plié A t = {i.e. sur le filament non-plié 
Xi, . . . Xiooo); la plus grande valeur propre de la différent ielleF^ de est 
négative et comprise entre —5,2 et —3,5 (figure [H]). Elle atteint un minimum 
en Xqj5, c'est-à-dire à proximité du futur maximum de courbure Xggo- 

La seconde valeur propre est très éloignée de 1, et reste comprise entre 0,04 
et 0,07. Le premier vecteur propre correspond donc à la direction instable. 
Il est représenté en différents points du filament sur la figure [32J On peut 
le décomposer en trois parties : (1) < t < 1,4; (2) 1,4 < t < 1,6; (3) 
1,6 < t < 2. Il faut tenir compte de la normalisation du vecteur propre pour 

25. On a calculé ici les valeurs et vecteurs propres de la différentielle D de T^, et non 
ceux de D* D, qui permettent d'évaluer exactement les directions contractées ou dilatées, 
dans la mesure où < x, D*D{x) >= Ces deux résultats ont cependant de bonnes 

chances d'être très semblables. 
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Population mature au niveau du pli 

7| 1 1 1 1 




Ql 1 1 1 1 1 1 

-2-1 1 2 3 4 

t 



Figure 39 - Évolution en temps continu au niveau du pli. 



Vecteur tangent au filament non-plié, en x„„„ Vecteur tangent au filament plié, en T^Cx^^^) 




(a) en xego {t = 0) (b) en T^ix^go) {t = 2) 

Figure 40 - Vecteurs tangents au filament. 
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Plus grande valeur propre de la différentielle au niveau du pli (en x.) 




j 

Figure 41 - Différentielle de sur le filament non-plié : première valeur 
propre. 



interpréter correctement son évolution le long du pli. La première partie 
varie très peu le long du filament, ce qui confirme l'hypothèse de formation 
du pli à t ~ 1,4 que l'on a formulée d'après la figure l38k. La troisième partie 
est rectiligne et remonte en même temps que le pic de la deuxième partie. 
Elle traverse en Xqsq, à l'instant précis où le pic atteint son maximum. Le 
pliage semble ainsi correspondre à une forte expansion localisée en t = 1,5, 
accompagnée d'une absence d'expansion (ou plutôt une contraction puisque 
les autres vecteurs propres ont de très petites valeurs propres) sur l'intervalle 
[1,6; 2]. Ces deux phénomènes sont difficiles à distinguer dans la mesure où 
la normalisation du vecteur propre entraîne une compensation entre la 
hauteur du pic et celle de la partie rectiligne. On reconnaît tout de même une 
forme de pliage en nous limitant à la partie (3) du vecteur propre, l'expansion 
ayant lieu dans des sens opposés de part et d'autre du pli. 

On remarque également que la direction instable à t = (figure H2|) est 
semblable au vecteur tangent ht = 2 (figure HUb). Cela signifie que le filament 
est suffisamment transverse à la variété instable en Xqqq. 

Les quelques irrégularités observées sur les courbes des figure HT] et HSl 
semblent confirmer l'hypothèse avancée lors de l'étude de la courbure du pli : 
la ligne brisée considérée doit «sauter» d'un filament à l'autre, deux filaments 
côte à côte n'ayant pas précisément la même courbure. 
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Premier vecteur propre de la différentielle de en x=x. 
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Figure 42 - Différentielle de sur le filament non-plié : premier vecteur 
propre. 

Filament plié Àt = 2 {i.e. sur le filament plié T^(xi), . . . T^(a;iooo)), les 
deux plus grandes valeurs propres de la différentielle de sont représentées 
figure uni Elle est légèrement supérieure à 2 au voisinage de T^(x69o), mais il 
n'y pas d'expansion tout au long du filament. Celle-ci est en effet de module 
inférieur à 1 entre T'^{xi) et T^(x444) où elle traverse le cercle unité, puis 
croît lorsqu'on se déplace le long du filament plié vers Xiooo où elle dépasse 
4. Il n'y a donc pas uniforme hyperbolicité sur l'attracteur. 

Dans le domaine où les valeurs propres sont toutes de module strictement 
inférieur à 1, la courbe Xmax{j) n'est pas continue. En traçant aussi la se- 
conde valeur propre, on constate qu'il s'agit d'un échange dans l'ordre des 
deux premiers vecteurs propres. Tant que la première valeur propre n'est pas 
suffisamment grande devant la seconde, le premier vecteur propre ne suffit 
pas à décrire la différentielle de T^. Il est cependant intéressant de le considé- 
rer autour de T^(a;69o) (figure Hl]). Il a alors le même aspect qu'à la figure H2| 
avec un pic autour de t = 1,6. 

Compte-tenu du vecteur tangent au temps t = A (figure Hïïk). le filament 
plié est faiblement transverse à la variété stable en T^(a;69o)- Après le pliage, 
il semble donc que la principale action de soit de renforcer le pli, par une 
forte contraction, l'expansion étant très faible dans la direction tangente au 
filament en T^(x69o)- H y a par ailleurs assez peu de différences entre les diffé- 
rentes directions instables représentées figure SU Le pli étant déjà formé, il y 
a expansion dans la direction tangente au filament loin de T^(a;69o)- En effet. 
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Valeurs propres de la différentielle au niveau du pli (en "r^(x.)) 



CO 2 
E 



Première valeur propre 
Deuxième valeur propre 



100 200 300 400 500 



700 800 900 1000 



Figure 43 - Différentielle de sur le filament plié : premières valeurs 
propres. 



Premier vecteur propre de la différentielle de en x=T^(x.) 
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Figure 44 - Différentielle de sur le filament plié : premier vecteur propre. 
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Vecteur tangent au filament plié, en T^(x ) Vecteur tangent au filarrent plié, en T*(x ) 




(a) en T^(x69o) (b) en T'^(x9oo) 

Figure 45 - Vecteurs tangents au filament après pliage (t = 4). 



le vecteur tangent en T^(a;9oo) (figure H5b) coïncide avec la direction instable 
au temps t = 2. C'est également le cas du vecteur tangent en T^(x4oo), qui 
est précisément opposé au vecteur tangent en T^(x9oo)- 

Visualisation animée de l'attracteur Une dernière façon de comprendre 
la formation du pli (et la dynamique générale de l'attracteur) est d'utiliser 
une animation^ représentant la projection dans de l'attracteur avec ori- 
gine des temps St pour des valeurs successives de St (voir section l4.2.ip . 

On constate pour —0,6 < ô < —0,35 qu'une portion jusque là quasiment 
rectiligne se plie littéralement au cours d'un rapide déplacement dans de la 
zone «inférieure» (deux coordonnées proches de 1, la troisième grande) vers la 
zone «supérieure» (une coordonnée proche de 1, les deux autres grandes). Le 
pli est inexistant pour ôt = —0,6 (figure H^). quasiment formé pour 6t = —0,5 
mais pas encore placé en x^go (figure H6b) et l'est totalement à ôt = —0,35 
{i.e. t = 1,65 avec les conventions de cette section : voir figure H6b) . 

La suite de la déformation, pour ôt > —0,4, ne fait qu'accentuer ce pli en 
étirant l'attracteur dans une direction et en le contractant dans les autres (au 
voisinage du pli). Cette contraction se ressent plus particulièrement autour 
de ôt = 0,5 (z.e. t ^ 2,5 : figure SSl). 

26. Il s'agit du fichier film_delta.avi. 
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N{2l+2=0.6) 1 N(2t+1-0.6) 

(a) ôt = -0,6 



N(2t+2-0.5) 



1 1 

(b) ôt = 



N(2t+1-0.5) 

-0,5 
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Figure 47 - Spectre de la différentielle de en les 80 points de la carte. 



4.2.9 Non-hyperbolicité de l'attracteur 

La figure HT] représente le spectre de la différentielle de en chacun des 
80 points de la carte (en mettant la structure canonique sur chaque espace 
tangent). Cela nous permet de caractériser l'action de sur les différentes 
zones de l'attracteur. 

La figure SHl permet de situer dans M.^ les zones d'expansion (dans la ou 
les direction(s) instable(s)) et de «non-hyperbolicité» 0. 

L'expansion est assez forte dans l'ensemble, à l'exception de trois régions : 
le ph (25, 28, 22, 27), le creux (61, 62) et le triangle (78; 65; 75, 72, 73). 
Ces trois zones semblent ainsi jouer un rôle particulier dans la dynamique. 
Nous avons déjà commencé l'étude du pli, il faudra également s'intéresser 
aux deux autres zones. 

27. non-hyperbolicité avec la structure sur chaque espace tangent. Il est possible 
qu'une autre structure rende l'attracteur hyperbolique. 
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ExpEinsion et hyperbolicité (e^ = 0.17609) 




Figure 48 - Spectre de la différentielle selon les zones de l'attracteur. 
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Le «triangle» possède sans aucun doute une orbite de période 3 pour T^. 
Il faudrait évaluer plus précisément l'hyperbolicité de cette orbite, qui joue 
un grand rôle dans la dynamique globale. 

Le «creux» semble quant à lui être une autre zone de pli dans l'attracteur. 
Il faudrait déterminer si c'est effectivement le cas, et si ce pli est ou non 
distinct du pli que nous avons déjà mis au jour. 

Enfin, deux points présentent à la fois une forte expansion et un défaut 
d'hyperbolicité (c'est-à-dire une seconde valeur propre proche de 1) : 42 et 
70. Leurs images par étant proche du point 62, il est possible que leur 
particularité soit simplement liée au «creux». 

4.2.10 Stabilité de la structure dans l'espace des paramètres 

Après avoir détaillé certains éléments de la structure d'un seul attracteur, 
celui que nous avons observé pour {Aq, p, 7) = (0,15; 0,30; 8,25), on peut s'in- 
terroger sur la persistance de cette structure quand on fait varier légèrement 
les paramètres. Nous pouvons d'ores-et-déjà esquisser une réponse, visuelle, 
à l'aide d'une animati;nB représentant l'attracteur (A; 0,30; 8,25) dans 
quand Aq varie (voir aussi la section I4.1.4[ consacrée à cette exploration, 
notamment le diagramme de bifurcation [T^ . 

Pour 0,135 < Aq < 0,160, l'attracteur grandit continûment en deux par- 
ties, à partir d'une forme très simple (uniquement la région de l'équilibre et 
du pli, i.e. zones 2 et 3). Pour Aq ^ 0,139, les autres régions apparaissent 
brusquement, dans une configuration assez similaire à = 0,15. La crois- 
sance de l'attracteur se fait par l'allongement de filaments s'achevant pour 
une pointe, un peu de la même façon que lorsqu'on a observé comment la 
variété instable s'étend à l'intérieur de l'attracteur (voir section 14.2. 7p . L'es- 
sentiel de la structure est conservé, même si la complexité de la dynamique 
va apparemment en s'accroissant. 

En Aq = 0,16, les deux composantes connexes se rejoignent pour n'en 
former qu'une, qui continue à grandir de la même façon. Mises à part quelques 
interruptions (une ou deux fenêtres de périodicité), l'attracteur reste intact 
(tout en se complexifiant au fur-et-à-mesure de sa croissance) jusqu'à Aq ^ 
0,19. Il disparaît alors par brusques paliers pour laisser place à une grande 
fenêtre de périodicité, puis un cycle se refermant sur un point fixe. On ne le 
retrouve plus trace ensuite de cet attracteur. 

Un zoom sur la région du pli (zone 3) nous a permis de déterminer si 
la pointe où se situent les points 29 et 30 continue à grandir en se pliant, 

28. Voir les fichiers f ilm_A0_100__200_300_10_50_08250_l_2coul_0 . 1_0 . 3 . avi pour 

une vue d'ensemble, et f ilm_zoomA0_100 200_300_10_50_08250_l.avi pour un zoom 

sur la région pliée (zone 3). 
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parallèlement aux filaments des points 21 à 28. Contrairement à ce que l'on 
aurait pu penser, ce n'est pas le cas : la double pointe 29-30 reste d'un 
seul côté du pli sans le franchir. Il n'y a donc pas de complexification de la 
structure de cette manière-là. 

5 Perspectives 

Il reste beaucoup de travail à faire pour comprendre ce modèle. A la lu- 
mière des simulations dont nous venons d'analyser les résultats, nous pouvons 
cependant déjà esquisser quelques pistes afin de poursuivre ce travail. 

5.1 Simulations numériques 

Tout d'abord, plusieurs aspects restent à étudier numériquement afin 
de comprendre au mieux le comportement du modèle. Certaines questions 
étaient déjà évoqués dans | YB98] . d'autres sont venues à la lueur des nou- 
velles simulations numériques. 

Explorations 

- Lors de l'augmentation de Aq dans le diagramme (Aq, 0.30, 8.25), les 
deux filaments se rejoignent-ils préciséments ou bien se rapprochent-ils 
sufîisamments pour induire le mélange ? 

- Dans le cas d'une exploration particulière, chercher le plus grand nombre 
possible d'attracteurs, notamment en «suivant» automatiquement cha- 
cun jusqu'à la perte de stabilité, et en utilisant plusieurs conditions 
initiales. Ce travail pourrait déjà être fait dans le cas des explorations 
déjà faites, notamment pour éclaircir la question des discontinuités du 
diagramme. 

- Exploration plus exhaustive de l'espace des paramètres. 

- Tester d'autres transitions de l'équilibre stable vers la dynamique chao- 
tique. 

- Déterminer des régions de l'espace des paramètres (en 2 ou 3 dimen- 
sions) où se produisent les divers comportements observés (orbite pé- 
riodique attractive, cycle attractif, attracteur de type Hénon et autres 
attracteurs étranges), et les bifurcations qui se produisent à l'inter- 
face entre ces différentes zones. Pour chaque comportement observé, on 
pourrait chercher à caractériser un peu plus précisément la dynamique 
dans un cas particulier au moins (notamment pour les cycles et les 
attracteurs de type Hénon). 
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- Déterminer par un calcul numérique la nature des bifurcations obser- 
vées sur le diagramme, en calculant les valeurs propres de la diffé- 
rentielle (soit en linéarisant l'équation, soit expérimentalement). Ceci 
serait particulièrement utile au niveau des discontinuités apparentes du 
diagramme. 

Étude de l'attracteur étrange 

- Géométrie de l'attracteur : 

- Décrire géométriquement l'attracteur avec plus de précision, et no- 
tamment comment se séparent les différents filaments (en se limi- 
tant aux plus gros d'entre eux). Ceci devrait permettre d'élaborer 
un modèle simple pouvant générer un tel attracteur. Il faudrait en 
particulier caractériser les «embranchements». 

- Y a-t-il des zones plus «denses» que d'autres (au sens de la mesure 
de Lebesgue et non de la mesure physique) ? 

- Dynamique sur l'attracteur : il reste beaucoup de travail pour com- 
prendre précisément l'intéraction entre cette dynamique en trois di- 
mensions et la dynamique en temps continu. 

- Dimension fractale : est-elle la même dans toutes les régions de l'at- 
tracteur ? 

- Sensibilité aux conditions initiales : calculer les exposants de Lyapunov. 

- Équilibre et orbites périodiques : 

- S'assurer de la correction de l'estimation de la seconde plus grande 
valeur propre de la différentielle à l'équilibre, éventuellement en li- 
néarisant directement l'équation. 

- Localiser quelques orbites périodiques de faible période (en particu- 
lier l'orbite de période 3 pour T^, située dans la zone 8 et l'orbite 
de période 2 dans les zones 4-5) et y effectuer la même étude que 
pour l'équilibre, notamment en évaluant les valeurs propres de la 
différentielle de T™ et en traçant la variété instable. 

- Y a-t-il d'autres équilibres, associés ou non à ces orbites, situés en- 
dehors de l'attracteur? 

- Le tracé que nous avons fait représente l'ensemble a;-limite : où sont 
situés les points périodiques? sont-ils denses? 

- Hyperbolicité : utiliser les orbites périodiques présentes sur l'attracteur 
pour la tester avec plus de généralité. 

- Ph(s) : 

- Peut-on éviter les discontinuités de la courbure en augmentant encore 
la précision de localisation du pli ? 

- Y a-t-il un lien entre pics de courbure et défaut d'injectivité de la 
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projection ? 

- Localiser les autres plis s'il y en a (notamment en étudiant le «creux» 
et sa préimage, qui se sont distingués dans l'étude de l'hyperbolicité). 

- Mesure physique : Evaluer la mesure physique sur l'attracteur : quelles 
zones sont plus chargées, quelles zones le sont moins ? 

- Essayer de «suivre» précisément la structure mise en évidence sur cet 
attracteur lorsque l'on fait varier légèrement les paramètres (l'équilibre, 
le pli, les valeurs propres de la différentielle dans ces régions, etc.). 
Jusqu'où peut-on la suivre ? 

D'autres attracteurs étranges ? Le cas (0,18; 0,30; 8,25) (figure fT2|) semble 
plus complexe mais peut-être encore plus intéressant que l'attracteur que 
nous avons étudié. Parmi les explorations déjà effectuées (ou à venir), on 
pourrait chercher à approfondir l'étude de la dynamique de certains objets 
remarquables, à la lumière du travail déjà effectué. 

5.2 Conjectures 

Au vu des résultats numériques, on peut énoncer quelques conjectures, 
en vue de rendre rigoureuses les observations qualitatives que nous venons 
de faire, et dont nous sommes à peu près sûrs. 

La dynamique est chaotique pour certaines valeurs des paramètres. L'at- 
tracteur est étrange. La dynamique chaotique est persistante, au voisinage 
de cet attracteur étrange. 

La dimension fractale de l'attracteur étrange est comprise strictement 
entre 1 et 1,5. 

L'attracteur est «quasi-hyperbolique» (cf. propriétés d'un système hyper- 
bolique et de l'attracteur de Hénon en annexe [C.3I) . 

On peut le décomposer en un nombre fini de parties dans lesquelles 
est transitive («décomposition spectrale»). 

Il y a une cascade sous-harmonique lorsqu'on se déplace dans l'espace 
des paramètres à partir de 7 petit, p proche de ou 1, petit, pour se 
diriger vers des valeurs plus grandes de 7 ou de p, ou bien vers des valeurs 
intermédiaires de p. 

5.3 Questions biologiques 

Il y a de nombreuses manières de complexifier le modèle pour le rendre 
plus réaliste, par exemple en ne supposant plus que la survie est indépendante 
de l'âge. On peut également essayer d'introduire des effets maternels. 
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6 Conclusion 



Revenons tout d'abord au problème purement biologique que nous nous 
sommes posé initialement. Il s'est avéré qu'un modèle simple combinant ma- 
turation et des saisons régulières peut engendrer des comportements chao- 
tiques extrêmement complexes, pour des valeurs assez raisonnables des para- 
mètres. La principale condition, et qui se trouve remplie pour ce qui concerne 
Microtus epiroticus, est une très forte fécondité. 

Il est clair que ce modèle pourrait difficilement être prédictif, tant il a été 
simplifié sans se soucier finement du cycle de vie des populations qui nous 
intéressent. En revanche, nous pouvons d'ores-et-déjà tirer des conclusions 
qualitatives, la plus importante étant que l'imprévisibilité des effectifs futurs 
à long terme peut avoir lieu dans un environnement stable, où tous les hivers 
sont strictement identiques. 

De plus, de nombreuses difficultés techniques soulevées lors de la mise en 
œuvre de simulations sur ce premier modèle seront très utiles pour des calculs 
numériques sur des modèles plus complexes dérivés ou non de celui-ci. 

D'un point de vue mathématique, nous avons eu un aperçu de la richesse 

des comportements que peut générer un modèle somme toute assez élémen- 
taire. En nous attardant sur l'un de ces «attracteurs étranges», nous avons 
pu évaluer la complexité d'un seul de ces comportements, en mettant en évi- 
dence des phénomènes très mal compris, voire jamais abordés encore. Nous 
avons finalement posé beaucoup plus de questions que nous n'avons apporté 
de réponses. De nombreuses simulations restent ainsi encore à faire, soit pour 
confirmer une hypothèse expliquant les résultats obtenus, soit pour nous ai- 
der à en formuler au sujet des questions pour lesquelles nous n'arrivons même 
pas à esquisser une réponse. 

Le travail théorique restant à faire est lui aussi immense. Avec quelques 
outils élémentaires, nous avons pu définir l'attracteur global du système, 
mais nous n'avons aucune information à son sujet. Les expérimentations nu- 
mériques nous ont permis de formuler quelques conjectures, mais ont surtout 
pour objectif de nous indiquer par quelles voies il serait possible de s'attaquer 
à la résolution de celles-ci. Au vu des quelques figures que nous avons pu tra- 
cer, il semble assez probable que les «attracteurs étranges» — s'il s'avèrent en 
être effectivement — que nous avons observés ont certainement un grand in- 
térêt dans l'étude théorique des systèmes dynamiques non-uniformémement 
hyperboliques. 
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A Simulations : travail préliminaire 



A.l Discrétisation du modèle 



On veut passer du modèle continu décrit par l'équation (12. ip à un modèle 
discret, qui pourra être simulé numériquement. 

Fixons un entier p > 0, ce sera le nombre de classes d'âges considérées par 
année. On note le nombre de naissances qui ont lieu dans l'intervalle de 



temps 



i—l . i_ 
P ' P 



Ni l'effectif mature moyen au cours de ce même intervalle, Cj le 



facteur saisonnier égal à la moyenner I dans cet intervalle de mp ^ (définie par 
(12.8p . i étant considéré mod. p). Soit Si la proportion des individus matures et 
vivants parmi ceux qui sont nés i pas de temps auparavant, et m la fonction 
de fécondité définie par les relations (12.5p . On a alors les relations suivantes : 



^ m{Ni) xNiXei 
m = 



P 



2p 



(A.l) 



k=l 



La condition initiale ('T-î)i<j<2p étant donnée, ceci permet de calculer les 
fii pour tout i > 2p. 

Remarquons que le calcul de Sj n'est pas toujours évident. Pour les calculs, 
dans le cas où Aq est un multiple entier de 1/p {Aq = io/p), on a pris 

li>iQ. Lorsque Aq n'est pas adapté au pas de discrétisation. 



1 - 



pxAj 



on a pris pour Si la moyenne des Sj^q pour {i — 1) x k + 1 < j < i x k, on sj^g 
désigne le coefficient de survie calculé avec un pas q = px k (le plus souvent, 
k = 100). C'est le cas du diagramme [TU 

Pour des raisons pratiques de calcul, on a souvent utilisé le vecteur des 
naissances Ui au lieu des effectifs matures Ni, car cela évite de calculer deux 
fois Hi au cours de la simulation. Il est bien sûr très simple de passer des 
naissances aux effectifs matures, mais en perdant les deux premières années. 
Ce choix explique l'apparition des naissances dans les résultats, alors que 
celles-ci ne sont pas explicitement utilisées dans le modèle. 



A. 2 Choix des conditions initiales 

Pour choisir une condition initiale «au hasard», on a choisi de détermi- 
ner un vecteur de naissances aléatoire. Les naissances successives sont tirées 
suivant des lois uniformes indépendantes. Deux méthodes ont été employées. 

29. En pratique, on a pris la moyenne des valeurs aux extrémités de l'intervalle. 
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Condition initiale I : population mature de la première génération ; = 0.1 8 ; y = 8.25 




Naissances Première génération 



Figure 49 - Condition initiale (I). 



Pour la condition initiale (I), on a imposé qu'il n'y ait pas de naissances 
en hiver (avec p = 0,41), et que l'effectif mature soit égal à 20 à l'instant t = 
(avec Aq = 0,18, Ai =2). La figure Hïïl représente les naissances (—2 > t < 0) 
et la première génération de populations matures (0 < t < 2) qui en découle 
(avec Ao = 0,18, Ai = 2, 7 = 8,25). 

Pour la condition initiale (II), qui a été utilisée le plus souvent, on n'a 
pas tenu compte de l'hiver (le choix d'une valeur particulière de p n'étant ps 
justifié), et on a imposé un effectif mature égal à 1 à l'instant t = (cette 
valeur étant plus raisonnable au vu de la dynamique stationnaire du système). 
La figure [SÔ] représente dans les mêmes conditions que précédemment les 
naissances et la première génération correspondante. 

On pourrait bien sûr concevoir d'autres méthodes de choix aléatoire d'une 
condition initiale, évitant mieux les biais possibles, mais ce n'est pas très im- 
portant pour l'usage que nous en avons eu. Cela ne serait utile que dans le 
cadre d'une étude de la «taille» de bassins d'attractions de plusieurs attrac- 
teurs. 

A. 3 Choix du pas de discrétisation 

Le paramètre p (nombre de classes d'âge par année, appelé an pour plus 
de clarté) est décisif pour les simulations. La complexité de l'algorithme est en 
effet proportionnelle au carré de ce pas. Le choix an = 100 est le résultat d'un 
compromis entre rapidité du calcul et précision, après quelques simulations 
test. Cette valeur est de plus raisonnable pour ce modèle : au vu de l'ordre 
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Condition initiale II : population mature de la première génération ; = 0.18 ; y = 8.25 




Naissances Première génération 



Figure 50 - Condition initiale (II). 



de grandeur des paramètres Aq et p, il ne semble pas utile d'être plus précis. 
De plus, la situation biologique de départ étant discrète, il est inutile de 
considérer une échelle de temps inférieure à trois jours. 

A. 4 Lissage des fonctions 

Le choix des fonctions m{N) et mp{t) étant assez arbitraire et peu réa- 
liste, le lissage des fonctions a surtout été fait à titre préventif, pour que 
d'éventuelles discontinuités ou irrégularités dans la dynamique ne soient pas 
dûes au manque de régularité des fonctions utilisées dans le modèle. Il a 
été motivé par l'observation d'angles sur certains attracteurs qui semblaient 
lisses par ailleurs, au cours de simulations préliminaires. 

Il semble en réalité que ces changements n'ont pas modifié qualitativement 
les propriétés globales du système dans l'espace des paramètres. En revanche, 
il est certain que pour des valeurs fixées des paramètres, la plus infime mo- 
dification des fonctions utilisées dans le modèle peut modifier entièrement le 
comportement observé. 



III 



B Simulations : traitement des données 



B.l Visualisation de l'attracteur en 3 dimensions 

On représente {N{t),N{t + l),N{t + 2) pour les valeurs entières de t G 
[^min; ^max]- Ces valeurs entières correspondent précisément à la fin de l'été 
dans le cas tête = 0, au milieu de l'automne dans le cas général. 

Intervalle de temps choisi Le plus souvent, on a choisi de se limiter à 
19 001 < t < 20 000. La valeur maximale est choisie pour que le régime 
transitoire soit largement dépassé, et elle fixe la durée du calcul : avec 100 
pas par an, un calcul sur 20 000 ans se fait en à peu près 2 minutes avec un 
ordinateur récentE^. Avec t > 19 001, on dispose de suffisamment de points 
pour distinguer clairement le type d'attracteur (orbite périodique, cycle, type 
Hénon, étrange, etc.) tout en évitant le régime transitoire. 

Ces deux valeurs ont été testées sur quelques exemples de valeurs des 
paramètres, en s'assurant que le régime transitoire est très largement passé. 
Lors des simulations suivantes, on a vérifiéEïl que c'était toujours le cas. 

Visualisation de la mesure physique Pour limiter la taille des figures 
représentant l'attracteur (0.15,0.30,8.25), nous avons représenté une partie 
seulement des points que nous avons calculés, de telle sorte qu'il y ait de 
l'ordre de 20 000 points sur chaque graphique. Certaines zones étant beau- 
coup plus denses que d'autres, cette limitation a été faite dans des proportions 
différentes suivant les huit zones de l'attracteur que nous avons définies. Il 
devrait donc y avoir environ deux fois plus de points dans les zones 2 à 5, sur 
tous les graphiques où figurent les 80 points de la carte. Cette transformation 
n'a ainsi pas été effectuée sur la figure [T6l qui permet donc de se faire une 
idée plus exacte de la mesure physique sur l'attracteur. 

B.2 Diagrammes de bifurcation 

Afin de comprendre le rôle que jouent les paramètres du modèle dans la 
dynamique, nous avons réalisé des animations représentant les graphiques en 
3 dimensions, l'un des paramètres variant au cours du temps. Pour retrans- 
crire une partie de ces informations sur un graphique bidimensionnel, nous 
avons dû tracer des diagrammes de bifurcation. Le principe est le suivant : on 
porte en abscisse l'un des paramètres du système et N{t) en ordonnée, pour 

30. PC avec un processeur Athlon 2,4 GHz et 768 Mo de mémoire vive. 

31. Cela se détecte facilement sur la représentation en trois dimensions. 
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t entier, 19 001 < t < 20 000. On s'est contenté de 1 000 valeurs de t, au lieu 
de 10 000, car la projection unidimensionnelle de l'attracteur ne permet pas 
de bien saisir sa géométrie. On distingue tout au plus les ensembles finis, les 
ensembles fractaux et les ensembles continus. Une précision accrue n'aurait 
rien apporté. 



B.3 Injectivité de la projection 



Pour évaluer l'injectivité de la projection vr : M — )■ M , nous avons 
cherché à évaluer sup^-^^/^pj ^^^^i"^^*^)^^"^!/^!^^'^^"^ au voisinage des 80 points de la 
«carte» de l'attracteur (section [4.2.1l) . En chaque point, nous avons déterminé 
les éléments de la boule B de rayon r = 0.1 dans L^(]R^°^) centrée en ce point 

2^201 (^o)- Nous avons alors calculé sup^^^^^B ''""'"xS^-^aftOluf °' " 

Nous avons indiqué une deuxième information sur les graphiques ainsi 
obtenus : le nombre de points dans chaque boule. En effet, certains points 
sont situés dans des zones beaucoup plus denses que d'autres (au sens de 
la mesure physique), et cela peut fausser les résultats obtenus puisque nous 
ne conservons que le maximum sur les points de la boule. Il est normal 
d'obtenir un résultat plus élevé si la boule contient plus de points, puisque 
l'on risque d'avoir des points très proches dans par «accident», au vu des 
approximations que nous sommes obligés de faire. 

Nous avons choisi les normalisations suivantes pour les normes dans 
M^, p < oo : Il (xi, . . . , a;Ar)||^p = -j^ x^- La norme L°° est simplement le 
sup des coordonnées. 

Notons que lorsque ôt est pris non-nul, le vecteur X2oi{t) est décalé d'au- 
tant. Les boules considérées ne sont donc pas le mêmes lorsque ôt est différent. 
Ceci permet de mieux comparer les différentes valeurs de ôt, et de ne pas pri- 
vilégier ôt = 0. 

Pour résumer les nombreux graphiques obtenus, nous avons choisi d'ex- 
traire des 80 valeurs de sup deux données : le maximum et la valeur médiane. 
Elles nous permettent d'avoir une assez bonne idée de la qualité de la pro- 
jection. 

Enfin, lorsque nous avons considéré d'autres visualisations, par exemple 
avec une échelle logarithmique, nous avons procédé aux mêmes opérations, en 
remplaçant X3 par gixs), X201 étant inchangé. Il y a alors un inévitable chan- 
gement d'échelle homothétique : par exemple, avec g{x) = Ax et A très grand, 
le résultat serait artificiellement bon. Il faut donc rapporter les quantités cal- 
culées au diamètre de la projection de l'attracteur. Dans le cas de l'échelle 
logarithmique, cela ne fait que confirmer notre conclusion en sa défaveur. 
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B.4 Géométrie 



La représentation simplifiée de la figure [22] est simplement un extrait de 
la variété instable, n = 18 (voir annexe [Bl9|) . Le découpage a été fait morceau 
par morceau, et le choix des régions est purement visuel. 

B.5 Décomposition spectrale, mélange 

Pour déterminer s'il y a ou non mélange topologique, nous avons utilisé 
plusieurs couleurs suivant la valeur de t modulo A^, pour quelques valeurs de 

entre 2 et 10. Lorsque des zones bien séparées se sont distinguées, nous 
avons conclu à la A^-périodicité de la dynamique. Dans le cas contraire, si 
nous avions l'impression d'un mélange des couleurs pour toutes les valeurs de 

(comme à la figure [12]), nous avons conclu au mélange topologique. Il ne 
s'agit donc que d'impressions visuelles, et non de vérifications rigoureuses. 

B.6 Calcul de la dimension fractale 

La dimension fractale est définie en annexe IC.9I 

Calculer la dimension fractale d'un ensemble dont on ne connaît qu'un 
nombre fini de points, avec une précision limitée, est loin d'être un problème 
facile. Nous avons donc dû faire de nombreuses approximations pour tenter 
d'obtenir une valeur approchée raisonnable. 

Attracteur (0,15; 0,30; 8,25) Nous avons considéré l'attracteur dans M^, et 
non dans R^°^, avec la projection naturelle, l'origine des temps étant prise 
ôt = 0. Il s'agit donc d'un ensemble K d'environ Aq = 100 000 points. 
Pour différentes valeurs de e, nous avons calculé le nombre N^{K) de cubes 
Ci,j,k = [^e; {i + l)e] x [je; (j + l)e] x [ke; {k + l)e] qui contiennent au moins un 
point de K. La figure [26] représente ainsi logiQ N^{K) en fonction de log^oe. 
En théorie, la dimension fractale est l'opposé de la pente limite en — oo 

de cette courbe (rigoureusement, on sait que 

pente limite doit donc être la même). Mais dés que e est assez petit, A^^ est 
constant égal à A'q = 100 000, car l'ensemble K est fini. Nous avons donc 
considéré la pente (obtenue par régression linéaire) en nous restreignant à 

< ^ = A"eo et e > y/eô. Le choix de ce domaine a été fait au vu des 
données et correspond à la zone où les points sont bien alignés. 

En raison de ces nombreuses approximations, il ne faut pas accorder trop 
d'importance à la valeur précise que nous avons obtenue, seul l'ordre de 
grandeur (entre 1 et 1,5) a de bonnes raisons d'être exact. 
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Diagramme (0.18,0.41,7) Pour réaliser la figure El nous avons dû calcu- 
ler la dimension fractale d'un grand nombre d'ensemble, via un processus 
automatisé. Au lieu de 100 000 points, nous avons dû nous contenter de 
Nq = 10 000 points. La précision de ces calculs est donc encore inférieure. 
Nous disposons probablement d'une légère sous-estimation des dimensions 
fractales de ces attracteurs. 

B.7 Sensibilité aux conditions initiales 

On choisit un point sur l'attracteur, considéré dans M^*^^, et l'on déter- 
mine l'ensemble des points de l'attracteur situés à une distance inférieure à 
r pour la norme L'^ (sur la figure [271 ^ = 0,04 et p = oo). Par point «de 
l'attracteur», on entend qu'il s'agisse d'un des 20 000 points d'une orbite 
calculée précédemment et dont on a les coordonnées dans M^*^^. Ces points 
sont donc répartis selon la mesure physique. 

Pour chacun de ces points, on trace N{t) — A^*(t) pour t G [—20; 20], 
où N*{t) correspond au point initialement choisi. L'intervalle de temps de 
référence est [0;2], le graphique de la figure [271 n'est donc pas exactement 
centré sur le présent. 

B.8 Point fixe 

Première localisation dans A la recherche d'éventuels points fixes, 
nous avons cherché à minimiser la distance (non- normalisée), dans R'^, 
entre x^lt) et X3{t + 2) = T'^{xs{t)). Sur les 50 000 points considérés (nous 
nous sommes limités à la «grande» composante), deux sont écartés de leurs 
images de moins 4.10"'^. Ces points étant très proches l'un de l'autre, nous 
avons considéré qu'ils sont à proximité d'un unique équilibre. La visualisation 
de ce point de son image en 3 dimensions ont confirmé cette impression dans 
la mesure où l'image du point «fixe» en est très proche et sur le même filament 
(figure [5l|) . 

En effet, d'autres points sur l'attracteur sont également à une petite dis- 
tance (inférieure à 2.10"^) de leur image par (à proximité du point 43), 
mais celle-ci est sur un filament bien distinct. Il n'y aurait donc qu'un seul 
point fixe sur l'attracteur. 

Localisation plus précise, dans R^°^ Afin de disposer d'une meilleure 
approximation de ce point fixe, aussi bien dans que dans R^°^, nous avons 
eu recours à une méthode un peu plus sophistiquée. 

Nous avons tout d'abord choisi deux points de part et d'autre de l'équi- 
libre présumé, avec une marge assez importante, et nous avons tracé le seg- 
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' '-k équilibre 

préimag|# . , \ ■ 



Figure 51 - Les «points fixes» et leurs préimages par T^. 



ment rejoignant ces deux points. Nous avons calculé les images successives de 
ce segment par T^, en les tronquant de telle sorte que l'on reste dans un même 
voisinage de l'équilibre présumé. Après un nombre suffisant d'itérations de ce 
processus, nous disposons d'une bonne approximation de la variété instable 
locale à l'équilibre. C'est en effet ce que montre le A-lemme sous certaines 
conditionso Le segment initial ayant été choisi transverse à la variété stable, 
et la troncature nous permettant déviter la région de pli, cela devrait être 
vrai dans notre cas pratique. 

Disposant de la variété instable locale, son image par T~^" doit se contrac- 
ter autour de l'équilibre pour n assez grand. Nous avons ainsi commencé par 
calculer T"^" de la variété instable, suffisamment tronquée pour éviter le pli, 
tout en gardant la mémoire du passé (les troncatures successives nous obli- 
geant à rajouter des points sur la ligne brisée). Pour une valeur suffisamment 
grande de n (de l'ordre de 10), nous avons obtenu, en appliquant T"^" à la 
variété instable locale obtenue, un morceau de variété instable très proche de 
léquilibre. 

En minimisant la distance entre les points de cette variété instable lo- 
cale et leur image (dans M^°^), on obtient ainsi une très bonne approximation 
de l'équilibre. L'erreur que nous avons obtenue est en effet environ égale à 

32. Le A-lemme ne s'applique sans doute pas directement dans notre cadre non hyperbo- 
lique, en dimension infinie, mais nous donne de bonnes raisons de penser que sa conclusion 
est au moins approximativement vraie. 
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10 ^, ce qui est presque 3 000 fois mieux que notre première approximationE§. 

B.9 Variété instable 

On part de la variété instable locale à l'équilibre. La variété instable 
globale est alors donnée par la définition lC.21[ Comme il s'agit d'un équilibre, 

11 s'agit de calculer les images successives par de la variété instable locale. 

Afin de garder une précision finale correcte, on a augmenté progressive- 
ment le nombre de points définissant la variété instable. On a rajouté des 
points intermédiaires entre deux sommets consécutifs de la ligne brisée dès 
que leur distance dans M.^ est inférieure à rj = 10~^. 

On a alors deux manières de visualiser la croissance de la variété instable 
dans l'attracteur. D'une part, les images successives T'^^{Wl^^{xeg)) pour 1 < 
n < 18. D'autre part, en ne considérant que n = 18, on peut faire grandir la 
variété instable en partant du voisinage de l'équilibre. L'ordre d'apparition 
des différentes zones est en principe le même. La première méthode est plus 
naturelle et plus simplement interprétable, mais la seconde donne plus de 
détails sur l'ordre dans lequel la variété instable se déploie dans l'attracteur. 

B.IO Calcul de la différentielle et de ses valeurs propres 

Etant donné un point x de R^°^, il est aisé d'estimer la différentielle de 
en ce point : on fixe e = 10~^ et on calcule (T^(x + ee^) — T^(a;))/e. Cela 
nous donne la dérivée partielle dans la direction Cj (de la base canonique 
de M^°^. La matrice des dérivées partielles nous donne une expression de la 
différentielle. 

Le calcul des valeurs propres et des vecteurs propres est alors réalisé à 
l'aide des fonctions intégrées de Matlab®. Cela signifie que les espaces tan- 
gents en X et T'^{x) sont égaux à l'espace euclidien ]R20i_ norme considérée 
est donc la norme L^(M^°^). Pour différentes valeurs de e (allant de 10"'^ 
à 10~^), les valeurs des cinq premières valeurs propres changent très peu : 
Al = -2,2942, Aa ^ 0,0433, A3 ^ -0,0283, A4,5 ^ 0,0214 ± 0,0028i Une in- 
certitude persiste cependant sur les valeurs propres A^, j > 2, dans la mesure 
où les algorithmes de calculs sont assez instables. Il serait utile de vérifier 
les résultats ci-dessus par un autre calcul de la différentielle, par exemple en 
linéarisant directement l'équation. 

Les vecteurs propres représentés sont normalisés pour que leur moyenne 
(L^) sur [0; 2] soit égale à 1. 

33. dans la mesure où la norme utilisée n'est pas normalisée 
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B.ll Hyperbolicité 



Pour le calcul de la différentielle et de son spectre, on a utilisé la méthode 
précédente avec e = 10"^. On a considéré qu'une valeur propre met en défaut 
l'hyperbolicité de l'attracteur lorsque (llog^ol-^)! < = logio(3/2)). Notons 
que les valeurs propres dépendent fortement de la structure des espaces tan- 
gents d'arrivée et de départ. Ce n'est pas parce que la structure canonique 
du fibré tangent ne rend pas l'attracteur hyperbolique que celui-ci ne l'est 
pas. 

Pour certains points, les deux premières valeurs propres sont égales (en 
module) : il s'agit des points où l'on a deux valeurs propres complexes conju- 
guées de module maximal. Dans les autres cas, il s'agit toujours d'une valeur 
propre réelle. 

B.12 Pli : Localisation 

La première étape dans la localisation du pli est arbitraire et un peu 
imprécise : on choisit deux points sur l'attracteur Xi et Xiooo (figure 15^). à 
peu près sur le même filament et tels que T'^{xi) et T^(xiooo) sont situés de 
part et d'autre du pli (figure [36b). 

On détermine ensuite un segment dans M^°^ joignant T~^^{xi) et T~^°(a;iooo); 
sur lequel on place 1 000 points (y compris les extrémités). En réappliquant 
T^" à ces 1 000 points, on définit une ligne brisée {xj)j=i ^xom '■ c'est le fila- 
ment non-plié de la figure [36k . qui correspond au temps t = 0. Le temps t 
correspond ainsi à la ligne brisée (T*(xj))j=i...iooo- 

Les détails concernant les calculs de courbure sont donnés à la section 
suivante. 

Les calculs des différentielles ont été effectués comme indiqué section [B. 101 
La différentielle a été calculée à deux instants distincts : avant pliage {t = 0) 
et à l'emplacement du pli (t = 2). 

Précisons enfin que lorsqu'on parle du filament à l'instant t = (ou de 
sa courbure, etc.), il s'agit des Xj. Dans le cadre d'une visualisation d'un xj 
particulier, dans M^°^, on note l'intervalle de temps [— 2; 0] : c'est l'instant 
final qui est pris en compte. De même, le filament à l'instant s est composé 
des T^{xj), et correspond à l'intervalle [s — 2; s] dans M."^^^. 

Deuxième tentative On a réalisé une deuxième tentative de localisation 
du pli, plus précise, devant en principe supprimer les irrégularités de courbure 
que nous avons constatées. Pour cela, on utilise la variété instable globale au 
point fixe calculée précédemment. Pour n = 13, une partie de f"'{W'^{Xeq)) 
est située au niveau du pli. On a alors déterminé sa préimage par /^^ et 
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augmenté le nombre sommets de la ligne brisée pour disposer d'un filament 
de plus de 1 000 points dans la zone de pliage (en recalculant son image par 



B.13 Vecteur tangent, courbure d'une ligne brisée 

Soient Xi, . . . , xjv les sommets successifs de cette ligne brisée. 



Vecteur tangent Le vecteur tangent Tj en un point Xj a été calculé tout 
simplement à l'aide de la formule suivante : 



Tj = ^— X 201. 



Tj est donc normalisé pour avoir une moyenne {L}) égale à 1 (c'est un élément 
de M^oi). 



Courbure En supposant les sommets de la ligne brisée assez proches, on 
peut évaluer la courbure de la courbe qu'elle approche à l'aide d'une version 
discrète des formules continues définissant la courbure : 



ds^ 



T, = 



\\Xj+i - Xj\\L2 
Xj+l - Xj 

ds.j 



dsj 



Les points Xj sont considérés dans et non dans R^^^. La norme 
utilisée pour estimer l'élément de longueur dsj n'est pas normalisé, mais cela 
ne change rien pour le calcul de Kj. On estime de cette façon la valeur absolue 
de la courbure, et non la courbure proprement dite. 
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C Rappels de systèmes dynamiques 



L'objectif de cette section est de rappeler les notions les plus importantes 
de théorie des systèmes dynamiques et que nous avons évoquées précédem- 
ment. Pour plus de détails, on se reportera par exemple à |KH9 5J ou [BP03] . 
Pour les aspects plus spécifiquement hyperboliques, on se reportera à |PT93] 
et |Yoc95] . 

Cl Un peu de vocabulaire 

Un système dynamique continu est la donnée d'un espace X et d'un 
groupe à un paramètre de transformation (c'est-à-dire une famille 

d'applications continues X ^ X telles que = o f^)- Un système 

dynamique discret est la donnée d'un espace topologique X {V espace des 
phases) et d'une application / continue X — )■ X. On peut se ramener du 
premier cas au second par le biais de l'application de retour de Poincaré 
Dans la suite, sauf indication contraire, on se placera toujours dans le cas 
discret. 

Définition Cl (Orbite). Si x est un point de X, V orbite (positive) de x est 
l'ensemble {f"{x)/n > 0}. Si / est bijective, V orbite de x est {f^{x)/n G Z}. 

La théorie des systèmes dynamiques s'intéresse particulièrement au com- 
portement des orbites. Il est souvent utile de considérer aussi des «pseudo- 
orbites» au sens de la définition suivante. 

Définition C.2 (5-pseudo-orbite). Soit 5 > 0. Une suite {Xn)nen (ou {xn)nei) 
est une 5-pseudo-orbite si Vi G N (ou Z), (i(/(xj), Xj+i) < 5. 

Par exemple, une orbite simulée numériquement est une (5-pseudo-orbite, 
puisque les calculs sont effectués avec une précision limitée. 

Définition C.3 (Partie invariante). Une partie A de X est dite invariant^^ 
par / si fiA) = A. 

Dans la suite, on supposera X métrique compact. 

La notion de conjugaison topologique est extrêmement importantes. Elle 
traduit l'idée que deux systèmes dynamiques sont topologiquement équiva- 
lents. 

Définition C.4 (Conjugaison topologique). Soit r > 0. Deux applications 
C' f : X X et g : Y ^ Y sont topologiquement conjuguées lorsqu'il existe 
un homéomorphisme h : X ^ Y tel que f = h^^ o g o h. 

34. attention, cette définition est parfois remplacée par f{A) C A. 



XII 



Lorsque h est un C"* difFéomorphisme (m < r), on parle de conjugaison 
lisse. Parfois, on peut seulement trouver h : X ^ Y continue surjective telle 
que h o f = g o h. On parle alors de semi-conjugaison. 

Nous pouvons désormais définir la stabilité structurelle d'un système dy- 
namique. 

Définition C.5 (Stabilité structurelle). Une application / est struc- 
turellement stable (1 < m < r < oo) s'il existe un voisinage U de f pour la 
topologie telle que toute application g E U est topologiquement conju- 
guée à /. 

Si de plus on peut choisir h = hg dans la conjugaison de f et g tel que 
hg et hg^ convergent uniformément vers l'identité lorsque g converge vers / 
pour la topologie C™, alors on dit que / est C" fortement structurellement 
stable. 

C.2 Récurrence 

Définition C.6 (Point périodique). Un point a; de X est dit périodique 
lorsqu'il existe n > 1 tel que /"(x) = x. On note Per{f) l'ensemble des 
points périodiques. 



/(Per(/)) = Perif) et /(Per(/)) = Per{f). 

Définition C.7 (Point récurrent). Un point a; de X est dit positivement 
récurrent (resp. négativement récurrent) si x est un point d'accumulation 
de la suite (/"(x))„>o (resp. (/~"(a;))n>o)- On note R^{f) l'ensemble des 
points positivement récurrents, R~{f) l'ensemble des points négativement 
récurrents et R{f) — R'^if) U R~{f) l'ensemble des points récurrents. 

Autrement dit, partant d'un point récurrent, on revient une infinité de 
fois dans son voisinage. On démontre : 

Per(/)ciî+(/)niî-(/)^0 
/(iî+(/)) = R^if) et /(/?-(/)) = R-{f) 

Définition C.8 (Point limite). Pour tout point x, on note uj{x) {ensemble 
uj-limite de x) l'ensemble des points d'accumulation de {f"'{x))n>o et a{x) 
{ensemble a-limite de x)l'ensemble des points d'accumulation de {f^{x))n<o- 
On définit alors V ensemble u -limite L^{f ) = Uxgx'^I^)' ensemble a-limite 
L~{f) = Uxex"(^) et Vensemble limite L{f) = L+{f) U L-{f). 
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Un point de uj{x) est un point dont l'orbite issue de x visite le voisinage 
une infinité de fois. On a les propriétés suivantes : 

R+{f)cL+{f) et R-{f)cL-{f) 
f{L^f))=L^f) et/(L-(/))=L-(/) 

Définition C.9 (Point errantQ). Un point x G X est dit errant s'il possède 
un voisinage V tel que f^{V) (IV = ^ pour tout n > 1. Sinon, on dit que x 
est non-errant, et on note fl{f) l'ensemble des points non-errants. 

Un point est donc non-errant lorsque tout voisinage se recoupe au moins 
une fois dans le futur. On peut démontrer les résultats suivants : 

TaïT) c n{f) et IFiJ) c n{f) 
f{nif)) = n{f) 

Définition C.IO (Point récurrent par chaîneQ. Un point est récurrent par 
chaîne si pour tout 5 > il existe une 5-pseudo-orbite périodique issue de x. 
On note C{f) l'ensemble des points récurrents par chaîne. 

Un point récurrent par chaîne est un point qui peut revenir exactement 
en lui-même en autorisant des «erreurs» d'amplitude aussi petites que l'on 
veut. 

/(C(/)) = C{f) et n{f) c C{f) 
C(/|a) = a si a g {P^),R{f),LU),C{f)} 

Proposition Cl. En résumé, on a : 

Perif) C R^if) C L±(/) C L{f) C n{f) C C{f) 

et chacune de ces inclusions peut être stricte. De plus, chacun de ces en- 
sembles (ainsi que leurs adhérences) est une partie invariante par f , au sens 
de la définition I C. 31 

Définition C.ll (Transitif). Un homéomorphisme / est transitif si pour 
tout ouvert non-vide U, UneN/"(^) dense dans X. 

Ceci équivaut à dire qu'il existe x E X dont l'orbite est dense {i.e. u{x) = 

X). 

35. en anglais : wandering 

36. en anglais : chain-recurrent 
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Définition C.12 (Mélange topologique). Un homéomorphisme / d'un es- 
pace métrique compact X est topologiquement mélangeant si pour tous U, V 
ouverts non-vides, il existe no G Z tel que Vn > uq, f^{U) fl 7^ 0. 

Si / est topologiquement mélangeant, alors est transitif pour tout 
entier 7^ 0. La réciproque est fausse (voir l'exemple de la rotation du 
cercle : IC.ip . Il n'y a en effet pas nécessairement de «mélange» de l'espace 
des phases sous l'action de la dynamique dans le cas d'un système transitif. Il 
y a ainsi transitivité de tout système restreint à une orbite périodique, mais 
jamais de mélange. 

Définition C.13 (Minimalité). Une partie fermée A C X est minimale 
pour / si elle est non- vide, invariante par / et si A ne contient pas de fermé 
non-vide invariant par / autre que A. 

Ceci équivaut à dire que l'orbite (positive) de tout point x G A est dense 
dans X. En particulier, tout point de A est récurrent. 

Exemple Cl. Le cercle 5^ est minimal pour la rotation Ra : 6 ^ 6 + a (mod. 
1) si a est irrationel. Ce système est transitif, mais pas topologiquement 
mélangeant. 

Exemple C.2. Le doublement de l'angle 6 ^ 26 (mod. 1) sur 5*^ est topolo- 
giquement mélangeant, donc transitif. 

Proposition C.2. Si X est métrique compact non-vide, f continue X ^ X , 
alors X contient un fermé minimal pour f . En particulier, R{f) 7^ 0. 

Les définitions suivantes précisent les notions intuitives d'attracteur et de 
bassin d'attraction. 

Définition C.14 (Attracteur). Une partie compacte A C X est un attrac- 
teur pour / s'il existe un voisinage V de A et un entier X G N tel que 

f^iV)cVetA = nneNny)- 

Définition C.15 (Bassin d'attraction). Soit A un attracteur. Le bassin d'at- 
traction de A, noté B{A), est l'ensemble des points x G X tels que uj{x) C A. 

C.3 Dynamique hyperbolique 

Le cas particulier de la dynamique hyperbolique est extrèment impor- 
tant, notamment parce qu'elle est présente dans la plupart des systèmes 
dynamiques. Le comportement hyperbolique est, comme nous allons le voir, 
le comportement «typique» d'un système dynamique. Nous parlerons ici de 
systèmes dynamiques uniformément hyperboliques. 
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C.3.1 Définitions 



Commençons par traiter le cas d'une application linéaire. 

Définition C.16 (Application linéaire hyperbolique). Une bijection linéaire 
T : E = R" M" est dite hyperbolique s'il existe une décomposition E — 
Es ® Eu en somme directe de sous-espaces T-invariants {i.e. TEg = Eg et 
TEu — Ey) tels que, en notant S — T\e^ et U — T\e„, il existe n > 1 tel que 
< 1 et ||C/-"|| < 1. 

Cette définition est indépendante du choix de la norme, et est équivalente 
à dire que T n'a pas de valeur propre de module 1. 

On dit que ||-|| est adaptée à T lorsque n = 1 convient dans la définition 
précédente et si \/Xs G Eg et G E^, on a + = max{||a;5||, On 
appelle alors constante d'hyperbolicité la constante 

ch{T) ^ max{\\S\\,\\U-^\\} < 1. 

Dans le reste de cette section, on considérera X — M une variété fisse, 
munie d'une distance d, U un ouvert de M et f : U — > M un difi^éomor- 
phisme sur son image. 

Définition C.17 (Point périodique hyperbolique). Un point périodique p de 

/, de période n, est hyperbolique si Dfp : TpM — )■ TpM est une application 
linéaire hyperbolique. Son orbite est appellée orbite périodique hyperbolique. 

Sans perte de généralité, on peut se limiter au cas d'un point fixe. Le 
théorème de Grobman-Hartman affirme alors qu'au voisinage d'un point fixe 
hyperbolique, un difféomorphisme est topologiquement conjugué à sa diffé- 
rentielle. 

Théorème C.3 (Grobman-Hartman). Soit Q un ouvert de R^, f : Q ^ 
un C^- difféomorphisme local, Xq un point fixe hyperbolique de f et T — 
Df{xo) la différentielle de f en xq. Alors il existe des voisinages ouverts U 
de dans R^ et V de xq dans Vt, et un homéomorphisme H : U ^ V tel 
que, pour tout x dans U avec T{x) dans U , on a 

HoT{x) ^ foH{x). 

Généralisons cette notion au cas d'un ensemble invariant quelconque. 

Définition C.18 (Ensemble hyperbolique). Une partie /-invariante A (Z U 

est hyperbolique si pour tout a; G A il existe une décomposition TM^ = 
E^ © E^ et des constantes A < 1 < et une métrique Riemannienne sur M 
vérifiant les propriétés suivantes : 
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- Vx e A, TUE^,) = et TUE^) = E^^^y 

- Vx e A, ||T/3;|B||| < A et < A*""^ (les normes étant induites 
par la métrique de M). 

On peut alors montrer que les sous-espaces E^ et E!^ dépendent continû- 
ment de X, ont des dimensions localement constantes, et sont uniformément 
transverses (il existe «o > tel que pour tous x G A, ^ e E^, rj G E^, l'angle 
entre ^ et est au moins cto). 

Définition C.19 (Anosov). Un difFéomorphisme f : M ^ M d'une 
variété compacte M est appelé un difféomorphisme Anosov si M est hyper- 
bolique pour /. 

L'ensemble des difFéomorphismes Anosov sur M est un ouvert de C^(M, M). 
Pour déterminer si un ensemble est hyperbolique, en général, on regarde 
s'il vérifie la condition de cône suivante. 

Proposition C.4 (Condition de cône). Soit A une partie f -invariante, U 
un voisinage de A, a > 1. On suppose qu'il existe en tout point x & U une 

décomposition TM^ = E], © E^. 

Supposons que Wx E U,\/v E T^M, en posant v = V1 + V2 et w = T^fiv) — 
wi + W2 (décomposition sur E]. et E^), on a : 

\\v2\\ > ll^ill =^ ||'"^2|| > ctllwill et \\w2\\ > û;||'y2|| 

\\w2\\ > ll'î^ill =^ Ikill ^ Q;||u»i|| et \\vi\\ > Q;||i'2||. 
Sous ces conditions, A est hyperbolique. 

Bien sûr, cette condition est sufiisante, mais pas nécessaire. 
C.3.2 Propriétés fondamentales 

Définition C.20 (Expansivité). Un homéomorphisme / : X — )■ X est ex- 
pansif s'il existe une constante Ôq > telle que pour tous x,y E X , x y, il 
existe n e Z tel que d{f^{x), f^{y)) > Ôq. 

Proposition C.5 (Expansivité). La restriction d'un difféomorphisme à un 
ensemble hyperbolique est expansive. 

Le lemme de pistage est fondamental pour justifier la validité de simula- 
tions numériques, où l'on fait des calculs avec une précision limitée sur un 
système dynamique chaotique. En revanche, il ne garantit pas que les orbites 
pistant les pseudo-orbites sont typiques. Ainsi, pour l'application f : x ^ 2x 
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(mod. 1), une orbite calculée par ordinateur s'achèvera toujours en 0, car la 
condition initiale est donnée par un nombre fini de bits. L'ordinateur calcule 
ainsi toujours une vraie orbite, mais systématiquement attirée par 0, ce qui 
n'est pas le comportement typique du système. 



Proposition C.6 (Lemme de pistager^l). Si A est un compact hyperbolique 
pour f , alors il existe un voisinage U{A) de A tel que pour tout ô > 0, il 
existe e > 0, V(x„)ngz e-pseudo- orbite contenue dans U{K), il existe x G X 
vérifiant Wn G Z, d{f'^{x),Xn) < à. 

Les orbites sous l'action d'une petite perturbation de / sont ainsi proches 
des «vraies» orbites de /, ce qui nous amène à la question de la stabilité 
structurelle. En fait, on peut même montrer qu'en un certain sens, la stabilité 
structurelle est équivalente à la notion d'hyperbolicité. 

Théorème C.7 (Stabilité structurelle). Si A est hyperbolique pour f : U 
M, alors pour tout voisinage V G U de A et tout S > 0, il existe e > 
tel que si f : U X et dci{f\y, /') < e, il existe un ensemble hyperbolique 
A' = /'(A') C V pour f et un homéomorphisme /i : A' — )■ A avec dco{Id, h) + 
dco{Id, h^^) < ô tel que le diagramme suivant commute : 

A' A' 



A A 

De plus, h est unique si ô est assez petit. 

En particulier, les difféomorphismes Anosov sont fortement structurelle- 
ment stables (voir définition \C.5\) . 



C.3.3 Variétés stables et instables 

La définition d'un ensemble hyperbolique dégage au voisinage de chaque 
point deux directions, l'une stable, l'autre instable, provenant de l'étude du 
système linéarisé. Intuitivement (le théorème de Grobman-Hartman fait déjà 
un pas dans ce sens) , le système non-linéarisé devrait présenter le même type 
de décomposition, au moins localement : cela conduit à définir les variétés 
stables et instables. Le théorème suivant justifie leur définition, dans un cadre 
local. 

37. shadowing lemma 
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Théorème C.8. Soit X hyperbolique pour un difféomorphisme f : V ^ 
M, avec des constantes \ < 1 < /i. Alors, pour tout point x & A, il existe 
W''^{x) et W^{x), images de disques par des C^-plongements, appelés variétés 
stables et instables locales en x, telles que 

1. T^W(x) = E'^ et T^Wix) = E^. 

2. f{W%x)) C W'U\x)) et f-^W^x)) C W{f-^{x)). 

3. pour tout S > 0, il existe C{ô) telle que pour tout n G 

Vy e W^ix), dinx),ny)) < C(5)(A + y) 
yy e W^ix), d{n{x),n{y)) < C{5){ii-5)-^d{x,y). 

4- il existe 13 > et une famille de voisinages contenant une boule 
autour de X G A de rayon (3 telle que 

w^'(^) = {y|r(y)eO^"(x),neN} 

iy"(x) = {y|r"(y)eO;-n(,),neN} 

Les variétés locales stables et instables ne sont pas uniques, mais l'inter- 
section de W({x) et W2{x) contient toujours un voisinage de x. On peut ainsi 
définir indépendamment du choix des variétés locales les variétés stables et 
instables globales. 

Définition C.21 (Variétés stables et instables globales). 

oo 

w'{x) = u mw^irix))) 

n=0 

oo 

W'^ix) = y r(l^"(/-"(a;))) 



n=0 



On a alors une caractérisation topologique des variétés stables et instables 
correspondant à la notion intuitive. 

Proposition C.9 (Variétés stables et instables). 

W^x) ^{ye U/dinx),ny)) > 0} 

n— >+oo 

W-{x) = {ye U/d{mx),n{y)) —-^ 0} 

n— >-+oo 

Ainsi, deux variétés stables (resp. instables) globales sont disjointes ou 
égales. 



XIX 



C.3.4 Produit local, ensemble localement maximal 

Lorsque des variétés stables et instables globales en un point s'inter- 
sectent, il est intéressant de considérer leurs points d'intersection. Ceci amène 
à définir le produit local, à l'aide de la proposition suivante (qui découle en 
partie du théorème IC.Sp . 

Proposition C.IO. Soit x E A. On note W^{x) et W^{x) les boules de 
rayon e dans W''^{x) et W^lx). Alors il existe e > tel que pour tous x G A, 
W^{x) n W^{y) contient au plus un point [x,y], et il existe ô > tel que si 
d{x, y) < à avec x,y E A, alors W^{x) fl W^{y) ^ 0. 

Définition C.22 (Produit local). On dit qu'un ensemble hyperbolique A a 
une structure de produit local si pour e > assez petit, le point [x, y\ défini 
par la proposition IC. lOl appartient à A. 

En fait, cette propriété est équivalente à la notion suivante (que l'on 
distinguera bien de la notion d'attracteur, définition IC. 141) . 

Définition C.23 (Localement maximal). Soit A hyperbolique pour f : U ~^ 
M. S'il existe un voisinage ouvert de A tel que A = Ay := Hnez /"^(^)î '^^ 
dit que A est localement maximal. 

Dans ce cas, on montre que les points périodiques de / sont denses dans 
l'ensemble des points non-errants de /|a. 

Théorème C.ll. Soit A un ensemble compact hyperbolique. Il a une struc- 
ture de produit local si et seulement si il est localement maximal. 

Une autre notion importante reliée au produit local est celle d'intersection 
homocline. 

Définition C.24 (Intersection homocline). Soit p un point fixe hyperbolique. 
Un point g 7^ p est homocline à p si g G W'^{p) fl W'^{p). Cette intersection 
est dite transverse homocline si les deux variétés stables et instables s'inter- 
sectent transversalement. 

Les exemples les plus importants sont le pendule (en temps continu) et 
le fer-à-cheval (en temps discret). En général, un système possédant une 
intersection homocline présente une dynamique très complexe. A ce sujet, on 
pourra consulter |PT93] . 
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C.3.5 Décomposition spectrale 

Il s'agit d'un résultat important qui permet d'étudier la récurrence des 
orbites dans le cas d'ensemble localement maximaux hyperboliques. 

Théorème C.12 (Décomposition spectrale). Soit A compact hyperbolique 
localement maximal pour un difféomorphisme f : U M. Alors il existe une 
famille finie des compacts invariants disjoints Ai, . . . , telle que ^^(/|a) = 
|J™^Aj. Les ensembles Ai) = Ai sont «irréductibles» au sens où /|a- est 

transitif. De plus, Ai = IJj=i ,i ^'^^^ fi-^i,j) ~ ^«j+i? f\A-i topologi- 
quement mélangeante. 

Si de plus A = C{f\\) {i.e. A est récurrent par chaîne) , alors le théorème 
suivant montre que A = fi(/|A), et donc la décomposition spectrale s'applique 
directement à A. 

Théorème C.13. Supposons que A = Per{f), R{f), L{f) ou C{f) est hy- 
perbolique. Alors cet ensemble est localement maximal et récurrent par chaîne. 
De plus, cet ensemble est égal à Per{f). 

Définition C.25. Si A est compact invariant par /, les ensembles stables et 
instables de K sont 

W'iA) = {ye X/ioiy) C A} 
W^{A) = {ye X/a{y) C A}. 

Cette définition est cohérente avec les notions précédentes de variétés 
stables et instables globales, comme le justifie la proposition suivante. 

Proposition C.14. Si A est un compact invariant hyperbolique localement 
maximal, alors W'{A) = UxeA^'l^^) W^"(A) = Ua;GA^"(^)- 

Définition C.26. On appelle ensemble basique l'un des Aj de la décompo- 
sition spectrale. Trois situations peuvent alors se présenter : 

- si W^*(Aj) est ouvert, on dit que Aj est un attracteur. 

- si W'^{Ai) est ouvert, on dit que Aj est un répulseui¥^ 

- si aucun de ces deux ensembles n'est ouvert, alors Aj est de type-selle. 

38. repeller en anglais 



XXI 



Tore T 



Image du tore = f(T) 




Figure 52 - Le solénoïde : T et /(T). 



C.3.6 Exemple de dynamique uniformément hyperbolique : le so- 
lénoïde 

Soit le tore T = {{9, z)/9 G M/Z, z E C, \z\ < 1} et l'application 
T > T 

f : 

(9, z) > {29 mod. 1, ^e^^'^ + ^z) 

Le tore et son image par / sont représentés figure [52j L'ensemble maximal 
invariant de / est A = flneN/"!-^)' hyperbolique^^, on l'appelle 

attracteur de Smale (figure [53k) . Localement, A est le produit d'un ensemble 
de Cantor dyadique (inclus dans le disque unité de M^) et d'une droite (fi- 
gure ISSb)- 

Les variétés stables sont les sections Cg^ = {{9,z) G A\9 = 9o}. Les 
variétés instables sont plus difficiles à décrire, et on peut montrer que chacune 
est dense dans A. 

Dans le cas du solénoïde, il est également possible de montrer l'existence 
d'une mesure physique (définition IC.33P qui donne la répartition statistique 
dans A des points de presque toutes les orbites. 



39. Les sous-espaces sont les plans 9 = cte, les EJ^ sont de dimension 1, tangents à 
A. 
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Solénoide : vue en coupe 




-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 



(a) A (b) coupe de A 

Figure 53 - Le solénoide. 



C.3.7 Dynamique non-uniformément hyperbolique 

Considérons l'exemple de l'attracteur de Hénon. Soit l'application 

^2 y M2 

H = Hb c : 

{x,y) > {x'^ + c — by,x) 

avec < 5 <^ 1 et c un peu plus grand que -2. 

Il existe un rectangle R tel que H{R) G R et H{R) ressemble à un arc 
de parabole «épaissi» (figure IMk ) : le rectangle est fortement pincé, étiré, 
et plié (dans le cas du solénoïde, il n'y avait pas de pli). L'attracteur est 
A = [j^ç^ H^[R), et possède en presque tout point une structure de Cantor 
X une droite (comme le solénoïde) (figure 15^). En revanche, il existe un 
ensemble de points (dense dans A) où ce n'est pas le cas : ce sont les «pointes». 
Ainsi, figure [5^ . on visualise une petite zone de l'attracteur, qui semblait 
rectiligne sur la vue d'ensemble, et on distingue une pointe, i.e. un filament 
qui ne se poursuit pas vers la gauche. 

On a alors : 

- une décomposition M."^ = (B E'^ en presque tous les points x G A, la 
contraction et l'expansion n'étant pas uniformes, et les espaces E^ et 
E^ ne dépendent pas continûment de x. 

- sensibilité aux conditions initiales. 

- le lemme de pistage reste valide pour presques toutes les pseudo-orbites 
(mais pas toutes). 
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(a) R et f{R) (b) A (c) zoom dans A 

Figure 54 - L'application de Hénon et son attracteur (6 = 0.02, c = —1.8). 



- pas de stabilité structurelle : avec des perturbations arbitrairement pe- 
tites de c, on peut obtenir une dynamique complètement différente. 
Cependant, pour presque toute perturbation, on a un attracteur du 
même type. 

- la mesure physique existe. 

C.4 Dynamique symbolique 

La dynamique symbolique fournit un exemple de système dynamique très 
important car il permet de modéliser la dynamique de très larges classes de 
systèmes. Nous la décrivons ici très brièvement. On trouvera une introduction 
plus complète dans |BP08] et |KH95] . 

Soit A un alphabet {i.e. un ensemble) fini, muni de la topologie discrète, 
et X = muni de la topologie produit, l'ensemble des mots sur A. C'est 
un espace métrique compact, muni de la distance 

d{uj,uj') = sup (2^1*1 fl^^^,^/) . 

On note cr : X — )■ X et on appelle décalage à (7anc/ie0 l'application définie par 
cr{ui) = Ui+i. Le système dynamique (X, cr) est appelé système de Bernoulli 
ou système symbolique. 

Soit i3 C un ensemble de règles de compatibilité. On appelle sous- 
décalage de type fini une partie de X 

S = {{uj.i)i(,z/{^h^i+i) e B,Wi e z}. 

E est invariante par a. 
40. shift en anglais 
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Ces notions permettent de coder de façon combinatoire la dynamique de 
certains systèmes, à l'aide de partitions de Markov. Une partition de Mar- 
kov d'un ensemble invariant compact hyperbolique localement maximal A 
est un recouvrement fini par des rectangles propres^ (Ra)a€A d'intérieurs 

o o 

disjoints et tels que si x G -Ra et f{x) G Rb, alors f{W^^{x)) C {f{x)) et 
W]^^{f{x)) C f{W]^^{x)). On définit alors un ensemble de transitions 

B = {{a,b)eA'/f{Ra)nRt^^ù} 

et S le sous-décalage de type fini associé. Le théorème suivant précise dans 
quelle mesure (S, a) «code» la dynamique de / sur A. 

Théorème C.15. Soient TZ = {Ra)a£A ^'^c partition de Markov de A pour 
f et (S, 0") le sous-décalage de type fini associé. 

1. Pour tout w G S, l'intersection flnez / "(-^^^n) ^^'^ réduite à un point 
7r(a;). 

2. L 'application tt : S — t- A esi continue, surjective, et le diagramme sui- 
vant commute : 



A A 

3. Pour toute mesure de probabilité a -invariante et ergodique fi, de support 
S, on a 

{{co G S/Card7r-i(7r(a;)) > 1}) = 0. 

Un ensemble basique Aj de la décomposition spectrale (Théorème IC.12P 
possède une partition de Markov de diamètre arbitrairement petit (voir |KH95] ). 



C.5 Chaos 

La notion de «chaos» en systèmes dynamiques, contrairement à sa signi- 
fication usuelle de désordre total, se réfère à une situation où les orbites ne 
convergent pas vers une orbite périodique ou quasi-périodique, et où l'évo- 
lution des orbites est imprévisible à un certain point, ou leur comportement 
est sensible aux conditions initiales. Les premiers exemples étudiés furent — 
entre autres — l'attracteur de Lorenz, l'application logistique et l'application 
de Hénon. 

41. Un rectangle R est une partie de A de diamètre plus petit que e/10 (e étant donné 

o 

par la proposition IC.10( ) et telle que [x, y] Q R si x,y ^ R. Il est propre si R = R. 
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Définition C.27 (Orbite chaotique). L'orbite de x, > 0}, est 

sensible (ou chaotique) s'il existe une constante C > telle que 



Une orbite asymptotique à une orbite périodique ou quasi-périodique n'est 
pas chaotique au sens où si f^^{x) et /"^(x) sont proches, alors /"^"'""(a;) et 
j"-2+"'(^x) restent proches pour tout n > 0. 

Une orbite sensible est également imprévisible dans la mesure où savoir 
qu'un point y de l'orbite est extrêmement proche de g G uj{x) n'est pas 
suffisant pour prédire le futur de y à une distance C près. 

Dans l'ensemble stable d'un attracteur hyperbolique non-trivial, de même 
que l'on a une forte sensibilité aux conditions initialeso, on peut montrer que 
l'ensemble des points ayant une orbite chaotique a une mesure de Lebesgue 
totale. 

Définition C.28 (Dynamique chaotique). Un système dynamique {X, f) est 
sensible (ou a une dynamique chaotique) lorsque l'ensemble des points ayant 
une orbite chaotique a une mesure de Lebesgue non-nulleF^. 

Cependant, le chaos ainsi défini ne doit pas être interprété comme une 
totale imprédictibilité. En effet, on observe numériquement, pour certains 
systèmes chaotiques, que pour toute condition initiale prise dans un certain 
ouvert, on obtient le même ensemble cj-limite. Ceci conduit à la notion d'at- 
tracteur étrange. 

Définition C.29 (Attracteur étrange). Une partie compacte A de X est un 
attracteur étrange s'il existe un ouvert U et N G U de mesure de Lebesgue 
nulle tel que Wx G U\N, u{x) = A et l'orbite de x est chaotique. 

Un exemple d'attracteur étrange est V attracteur de Hénon (section [C. 3. 7p . 
On appelle parfois également attracteur étrange un attracteur A tel que / a 
une dépendance sensible aux conditions initiales avec probabilité totale sur 
B{A) X B[A) (où B[A) est le bassin d'attraction de A : voir définition I Cl 5 1) . 

Une dernière notion importante est celle de dynamique chaotique persis- 
tante, qui traduit que de petites perturbations de / ont, avec une probabilité 
positive, une dynamique chaotique. Cette définition a un sens lorsque par 
exemple f = fa est paramétrée par a G M", car alors on dispose de la me- 
sure de Lebesgue sur l'espace des paramètre a. De façon plus restrictive, on 

42. Comme l'indique la propriété d'expansivité [CTSl 

43. Cette définition n'a de sens que lorsque X est une variété, pour que les ensembles 
de mesure de Lebesgue nulle soient définis. 



Wq G uj{x), Ve > 0, 3ni,n2,n > 0/ d{f'^^{x),q) < e, 

d{r%x),q) < e et d{r'+'' {x) , f'^+'^x)) > C. 
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peut demander la persistance d'une djTiamique chaotique dans un voisinage 
ouvert de /. 

Une notion que nous avons déjà introduite est étroitement reliée au chaos. 
Il s'agit de celle intersection homocline (définition IC.24p . Il y a équivalence 
entre l'existence d'une orbite chaotique (voir section \C.5\ définition IC. 2 7|) et 
l'existence d'une orbite homocline. En revanche, cela n'entraîne pas forcément 
que la dynamique est chaotique (définition IC.28p . 

Chaos et simulations numériques II est problématique de vouloir obser- 
ver ou même caractériser un comportement chaotique lors d'une simulation 
numérique. Comment en effet mettre en évidence un tel phénomène malgré la 
précision finie d'un ordinateur ? Celle-ci a plusieurs conséquences majeures. 

Tout d'abord, les erreurs d'arrondi font que l'on n'observe que des pseudo- 
orbites. Si le système étudié possède une propriété de pistage, comme c'est 
le cas avec les systèmes uniformément hyperboliques, on a de quoi être par- 
tiellement rassuré. Il reste cependant des cas (par exemple le doublement de 
l'angle) où les orbites qu'un ordinateur peut pister ne sont pas des orbites 
typiques du système. De même, lorsque les orbites calculées sont bornées, 
toutes les pseudo-orbites observées sont périodiques (même si la période est 
très longue), en raison du nombre fini de décimales que l'on peut calculer. Il 
faut donc fixer (arbitrairement) un seuil pour séparer orbites périodiques et 
non-périodiques. 

Un deuxième effet est que l'on ne peut observer que le comportement 
en temps fini. Comment alors être sûrs qu'il s'agit bien du comportement 
stationnaire, et non d'un régime transitoire très long? Il nous faut en effet 
fixer un seuil à partir duquel on observe la dynamique «à l'infini». Le choix 
de ce seuil est crucial pour éviter des erreurs, tout en limitant la durée des 
calculs. 

Enfin, lorsque l'on étudie un système dépendant de paramètres réels, il 
faut garder à l'esprit que l'on ne peut observer celui-ci que sur un ensemble 
de mesure nul, l'ensemble des rationnels. C'est tout l'intérêt de considérer la 
persistance de la dynamique dans un voisinage ouvert, Q étant dense dans 
R. Ce problème peut cependant se ramener à celui du lien entre pseudo- 
orbites et vraies orbites si la famille {fa)a£R dépend continûment de a pour 
la topologie de la convergence uniforme sur X, car alors une orbite sous fa+e 
est une pseudo-orbite sous fa si e est assez petit. 

C.6 Bifurcations 

Considérons une famille de systèmes dynamiques dépendant d'un ou plu- 
sieurs paramètres. Même si pour presque toutes les valeurs des paramètres, 
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le système a un comportement transverse (par exemple structurellement 
stable), il peut y avoir des valeurs particulières de ceux-ci où se produit une 
transition entre deux différents types d'orbites. De tels changements sont 
appelés bifurcations. Leur étude — qui est une branche à part entière de 
la théorie des systèmes dynamiques — est fondamentale pour comprendre 
les propriétés d'un système typique car les bifurcations montrent comment 
différents comportements transverses peuvent apparaître. 

Nous ne parlons ici que de quelques cas simples de bifurcations, en petite 
dimension. 11 en existe bien sûr beaucoup d'autres types. Nous nous limitons 
de plus à des bifurcations locales, c'est-à-dire pouvant être définies seulement 
au voisinage d'un point, par opposition aux bifurcations globales. Nous consi- 
dérons plus particulièrement le cas des bifurcations structurellement stables, 
définies de la manière suivante dans le cas de systèmes discrets. 

Définition C.30 (Bifurcation structurellement stable). Une famille {fr} de 
C°° difféomorphismes définis localement a une bifurcation structurellement 
stable à r = To si Z^-,, n'est pas localement structurellement stable et si pour 
toute famille {çr} de C°° difféomorphismes définis localement suffisamment 
C^-proche de {fr}, il existe une reparamétrisation 0(r) de {çr} et une famille 
continue {hr} d'homéomorphismes définis localement telle que 

94>{r) ^ fr^K 

partout où cela est défini. 

C.6.1 Diagramme de bifurcations 

Il existe un moyen simple de visualiser une bifurcation, appelé diagramme 
de bifurcation. On trace l'ensemble cu-limite L~^{f^) pour les différentes va- 
leurs du paramètre e, que l'on porte sur l'axe des abscisses. Un tel diagramme 
peut aisément être tracé numériquement, en prenant pour ensemble a;- limite 
les valeurs de /"(x) pour n «grand» et pour un ou plusieurs x choisis aléa- 
toirement. 

Il y a cependant une différence entre un diagramme obtenu par simula- 
tions et un diagramme théorique : les objets instables, ou de «petit» bassin 
d'attraction, n'apparaissent que dans le second cas. Il n'est ainsi pas for- 
cément simple de déterminer la nature d'une bifurcation en comparant son 
diagramme empirique avec les diagrammes théoriques des bifurcations clas- 
siques. 
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T <To 



T = To 



T >To 



Figure 55 - Bifurcation de la famille 



C.6.2 Cas discret, dimension 1 

En dimension 1, on peut classifier les bifurcations structurellement stables 
autour d'un point d'équilibre p. En effet, dans ce cas, la dérivée de fr^ en p 
doit valoir A = ±1. 

Commençons par le cas A = 1. La famille (/(+i).T)reiR) définie par 

Vx G M, = x + x^ + r (C.2) 

a une bifurcation structurellement stable en tq = 0, avec dérivée 1, et est 
caractéristique de cette situation. 



Proposition C.16. La bifurcation de la famille (1C.2|) en tq = est struc 



turellement stable, et toute bifurcation locale structurellement stable en di- 
mension 1 ayant lieu en un point fixe avec dérivée 1 est (topologiquement) 
équivalente (après reparamétrisation) à cette bifurcation. 

Ainsi, pour r < Tq, /(+i),t possède un point fixe stable —\/—t et un point 
fixe instable yj —t ; pour r = Tq, ces deux points fixes sont confondus, et 
l'équilibre qui en résulte est semi-stable; enfin, dès que r > Tq, /(+i),t n'a 
plus de point fixe (figure [55]) . Le diagramme de bifurcation correspondant est 
représenté à la figure [561 

Dans le cas où A = —1, le point fixe p est transverse et donc persistent. 
La valeur de la dérivée en p en supérieure à —1 pour r < Tq et inférieure à 
— 1 pour r > Tq, le point fixe restant isolé. Cela s'accompagne de la création 
d'une orbite stable de période 2, tandis que le point fixe devient instable. On 
parle de bifurcation par doublement de période., dont l'exemple typique est le 
suivant : 

/(-i),r(a;) = -TX + x^ (C.3) 
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Diagramme de bilurcation de fj^j ^ ^(xj = x^ + x + t 



équilibre instable 1 
équilibre stable | 




Figure 56 - Diagramme de bifurcation de la famille /(+i).r(3;), autour de 




T <To r = To r > To 

Figure 57 - Bifurcation subie par /(_i)^^(a;). 



au voisinage de a;o = 0, tq = 1. On montre alors une proposition similaire à 
la proposition IC.16| ce qui achève la classification dans le cas de la dimension 
1. Pour visualiser cette bifurcation, on peut tracer /(-i),t (figure [57]), mais 
aussi ^ (figure [58]) pour mieux comprendre les orbites de période 2. Le 
diagramme de cette bifurcation est représenté figure EHl 

C.6.3 Bifurcation selle-nœud 

En dimension supérieure, les bifurcations structurellement stables se pro- 
duisent lorsqu'une valeur propre vaut ±1 et les autres sont en-dehors du 
cercle unité. 

Un exemple classique, en dimension deux ou plus, est le suivant : deux 
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Diagramme de bifurcation de f, (x) = -t x + x' 



équilibre stable 

équilibre instable 
- période 2, stable 



0.5 0.6 0.7 



1 1.1 1.2 1.3 1.4 1.5 



Figure 59 - Diagramme de bifurcation de la famille ^-(x), autour de 
ro = 1. 
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i / J / / 
; / / / / 



1 \ \ \ \ 
^ \ \ \ \ 



f / 



fi=l> flQ 



1 < flQ 



Figure 60 - Bifurcation selle-noeud : espace des phases de part et d'autre 
de /io = 0. 



point fixes, un noeud (point fixe attractif) et une selle (attractif dans une 
direction, répulsif dans une autre) se rencontrent. Après bifurcation, il n'y 
a plus aucun point fixe (localement). Une telle bifurcation est appelée selle- 

Le système différentiel suivant donne un exemple de bifurcation selle- 
nœud : ^ 




X 



-17 = -y 



/i 



(C.4) 



Il ne s'agit en fait que d'une légère modification par rapport à la bifurcation 
\C2\ qui se produit sur la première coordonnée de ce système. La deuxième 
coordonnée est là pour que le point fixe instable devienne une selle (il ne 
peut pas y avoir de selle en dimension 1). Le diagramme de bifurcation est 
donc exactement le même que celui de la figure [56l L'espace des phases de 
part et d'autre de la bifurcation (yU = 0) est représenté figure [201 



C.6.4 Bifurcation de Hopf 

Etude d'un exemple dans M? Considérons l'exemple de la famille de 

systèmes dynamiques continus suivante : 




-\y + ex — ax(x^ + y^) 



= Xx + ey - ay{x^ + y^) 



(C.5) 



44. saddle-node en anglais 
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où A et a sont des constantes strictement positives. Pour tout e, (0, 0) est un 
équilibre du système, les valeurs propres de la dérivée en sont fi^ = iX + e 
et 717. L'équilibre est donc stable si e < et instable si e > 0. 

En coordonnées «polaires» (un peu modifiées), R = x'^ + y"^ et 6 = 
arctan^, (IC.Sp devient : 



- = 2iî(.-aiî) 

dO _ ^ 
'dt ~ 



(C.6) 



Ce système se résout explicitement (voir des exemples d'orbites figure lïïT] 
dans le cas a = 1, A = 2n), d'où 

- si e < 0, toutes les solutions convergent vers l'équilibre. 

- si e > 0, toutes les solutions (sauf la solution constante nulle) convergent 
vers l'orbite périodique 

a 

9 = \ 

La figure représente le diagramme de bifurcations de cette famille de 
systèmes dynamiques en eo = 0. 

Description du phénomène général La bifurcation décrite au para- 
graphe précédent est une bifurcation de Hopf. Plus généralement, considérons 
la famille à un paramètre d'équations différentielles dans 

Nous faisons l'hypothèse (Hq) : -fo(O) = et Do{Fq) n'a que des valeurs 
propres de partie réelle strictement négative, sauf deux qui sont imaginaires 
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Bifurcation de Hopf 



^— équilibre stable 
- - équilibre instable 
cycle stable 




Figure 62 - Diagramme d'une bifurcation de Hopf. 



pures et non-nulles : /io = et /ÏJJ = —iX, avec A > 0. Dans un voisinage de 
l'origine, le système peut se réécrire (après changement de variable), à des 
termes négligeables près : 

— — = — Axi — axQ[XQ + Xi) 

dx' , , 

— = Ax' 

y dt 

Nous supposerons désormais (Hi) : a > 0. 

Nous avons enfin besoin d'une dernière hypothèse relative à la dépendence 
en e. Dans un voisinage de 0, on peut suivre l'équilibre et les valeurs propres 
/ie, m proches de l'axe imaginaire. On suppose (H2) : ^3fî/ie > en e = 0. 

Sous les hypothèses (Hq), (Hi) et (H2), la dynamique de l'équation (1C.7P 
présente une bifurcation de Hopf au voisinage de l'origine en e = : 

- pour e < petit, il y a un équilibre stable. 

- pour e = 0, l'équilibre reste stable mais plus faiblement. 

- pour e > petit, l'équilibre est instable, mais une orbite périodique 
quasi-circulaire de rayon ^ A/e/a est stable. 

Cas des difféomorphismes Un phénomène semblable peut se produire 
pour des systèmes dynamiques discrets x ^ fe{x), x G M^. On fait les 
hjrpothèses suivantes : 
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1- /o(0) = 0, et les valeurs propres de -Dq/o ont toutes un module stricte- 
ment inférieur à 1 sauf deux, /io et JI^ pour lesquelles |/io| = 1- 

2. Pour k = 1,2,3,4, fi^ ^ 1, z.e. /iq ^ {±l,±i,±j}. 

3. (H'-l) et (H2) comme dans le paragraphe précédent. 

La dynamique pour e proche de est alors la même que dans le cas 
précédent. Un exemple de tel difféomorphisme est donné par 

f,{z) = X{l + e)z-az\z\\ z eC (C.9) 

avec |A| = l,A7^±l,a>0. L'équilibre est stable pour e < 0, faiblement 
stable pour e = cq = 0, instable pour e > et alors le cercle \z\ = est 
invariant et attire toutes les orbites proches de sauf l'équilibre lui-même. 

Remarquons également que si la dynamique sur la courbe invariante est 
proche d'une rotation, elle ne se comporte pas toujours comme une rotation. 
C'est le cas pour presque tous les paramétres, mais pas nécessairement pour 
tous. 

Cas des orbites périodiques On se ramène en fait aux difféomorphismes 
via V application de retour de Poincaré. En effet, soit l'équation différentielle 
^ = Fo(x) dans possédant une solution périodique xq- Considérons une 
section E transverse à l'orbite Xq en Xo(to)- Une condition initiale suffisam- 
ment proche de Xo(to) retourne sur E en un temps fini, ce qui définit (dans un 
voisinage de Xo(to)) un difféomorphisme /q de E. La même opération pouvant 
être faite pour une petite perturbation de Fq, cela définit une famille 
de difféomorphismes, comme dans le paragraphe précédent. 

C.6.5 Autres bifurcations 

Nous n'avons bien sûr pas abordé ici toutes les bifurcations possibles, 
même en nous limitant a priori à un cadre restreint. Un exemple particu- 
lièrement intéressant est celui de la bifurcation homocline, reliée à celle d'in- 
tersection homocline (définition IC.24P : deux intersections transverses homo- 
clines se rencontrent, forment une tangence à cet instant, puis disparaissent. 
Une référence à ce sujet est |PT93] . 

C.7 Dynamique des polynômes quadratiques 

La référence pour cette section est [Yoc99] . 
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On considère la famille d'applicationsri Pc : z t-^ z"^ + c pour z G C et 

ceC. 

Cette famille de systèmes dynamiques est l'une des plus simples qui, en 
dimension 1, peut générer un comportement chaotique. Son étude est de 
plus particulièrement intéressante car on y observe des phénomènes que l'on 
retrouve dans de nombreux autres cas. 



C.7.1 Ensembles de Julia et de Mandelbrot 

Il est intéressant de se placer dans C au lieu de M car on peut alors utiliser 
de nombreux résultats d'analyse complexe. Nous reviendrons ensuite au cas 
réel. 

Pour c G C, on considère V ensemble de Julia rempli (figure [63]) 
Kc = {ze C/P^{z) est borné} 
que l'on peut également écrire 

n>0 

en ayant posé R = {1 + a/1 + 4 |c|)/2. Ainsi, Kc est : 

- compact, 

- non-vide et il contient tous les points périodiques de P^, 

- totalement invariant, i.e. Pc{Kc) = Kc = P~^{Kc), 

- plein, i.e. C\Kc est connexe. 

Le bord Je = dKc est V ensemble de Julia. D'après un théorème montré 
indépendamment par Julia et Fatou, c'est aussi l'adhérence de l'ensemble des 
points périodiques répulsifs. 

Un théorème de Fatou (1919) montre que Kc si et seulement si 
Kc est connexe. Dans l'espace des paramètres, on considère V ensemble de 
Mandelbrot (figure 

M = {c E C/ Kc est connexe}. 

On montre que M = {c/ |P"(0)| < 2, Vn > 0} et donc M est compact. De 
plus, M est plein, symétrique par rapport à l'axe réel qu'il coupe suivant 
l'intervalle [—2,1/4]. Sur la figure [Mj on distingue des îlots disjoints de la 
grande composante de M. Un calcul plus poussé montrerait qu'ils lui sont en 
réalité reliés par des filaments extrêmement fins. 

45. Tout polynôme complexe de degré deux est conjugué par une application afHne à une 
application de cette forme. C'est en particulier le cas de la famille logistique x i— rx{l — x), 
bien connue en dynamique des populations. 
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Ensemble de Julia K ; c - =1 



Ensemble de Julia K ;c = 0.2 



Ensemble de Julia K ; c = 0.253 




t 


1 











-2 -1.5 -1 -0.5 6 1 1.5 2 -2 -1.5 -1 -0.5 0.5 1 1.5 2 



(a)c=-lGM (b)c = 0.2GM (c) c = 0.253 ^ M 

Figure 63 - Ensemble de Julia rempli K^. 



Ensemble de Mandelbrot 




Figure 64 - L'ensemble de Mandelbrot : approximation numérique. 
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Si c ^ M, Kc = Je est un ensemble de Cantor (figure [HSlc). La dynamique 
est de type décalage, et est hyperbolique. 

Si c G M, la présence d'orbites périodiques attractives {i.e. zq tel que 
P^{zo) = Zq et \{PJ!^y {zq)\ < 1) est déterminante pour la structure de Kc- 
Douady a montré (1982) que Pc a au plus une orbite périodique attractive. 
Lorsque c'est le cas, le bassin d'attraction W = {z/ lim„^oo d^P^iz), 0{zq)) = 
0} de l'orbite est l'intérieur de K^, et Pc\j^ est hyperbolique. Un tel c est alors 
dans l'intérieur de M, et la composante connexe de l'intérieur de M contenant 
c est appelée composante hyperbolique de M. 

Par exemple, l'ensemble des c tels que Pc possède un point fixe attractif 
est l'intérieur d'une cardioïde dite principale contenant 0. Pour un tel c, 
l'intérieur de Kc a une seule composante et Je est un quasi-cercle (figurel^b). 
La figure!^ a donne un exemple d'ensemble de Julia rempli lorsque c est dans 
une autre composante hyperbolique. 

L'intérieur de M est dense dans M et contient toutes les composantes hy- 
perboliques. La conjecture d'hyperbolicité dit que l'union des composantes 
hyperboliques est en fait exactement l'intérieur de M. Parmi les résultats 
partiels obtenus dans cette direction, on a montré que les composantes hy- 
perboliques de M rencontrent MflM = [—2, 1/4] suivant un ensemble dense. 

En revanche, cette intersection n'est pas de mesure totale, comme le 
montre le théorème de Jakobson |Yoc97] . 

C.7.2 Dynamique sur la droite réelle 

Supposons c & m n W = [—2, 1/4]. Nous venons de voir que pour un 
ensemble dense (mais pas de mesure totale) de valeurs de c, c est dans une 
composante hyperbolique et donc il existe une unique orbite périodique at- 
tractive, et son bassin d'attraction est l'intérieur de Kc- 

Partons de c = 1/4 et faisons diminuer c (voir le diagramme de bifur- 
cations, figure l^ïï]) . On a tout d'abord un point fixe attractif, puis une or- 
bite attractive de période 2 (après une bifurcation doublement de période 
en c = c^^^). Les bifurcations doublement de période se succèdent ainsi jus- 
qu'à atteindre c = c^°°^ où il n'y a plus d'orbite périodique attractive. Cette 
succession de bifurcations est appelée cascade sous-harmonique directe. 

Le point c = c^°°^ est appelé le point de Feigenbaum, où la dynamique peut 
encore être décrite assez simplement. Un intervalle / est stable, et possède 
deux sous-intervalles Iq et Ji disjoints tels que Pc{Iq) C /i et Pc{Ii) C /q. 
Dans chacun de ces intervalles /„, on trouve deux sous-intervalles 1^^ et 
Ia,i disjoints tels que P^^Iafi) C Ia,i et P^^Ia,!) C /a,o- On retrouve cette 
dynamique en faisant un double changement déchelle, spatial et temporel. 
L'ensemble limite a ainsi une structure d'ensemble de Cantor. Cependant, la 
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Diagramme de bifurcation de la famille quadratique 




Figure 65 - Dynamique de Pc pour c G [—2; 



dynamique de Pc sur cet ensemble n'est pas chaotique. 

Au delà du point de Feigenbaum, la dynamique devient chaotique. Il se 
produit alors un processus, miroir de la cascade directe, de regroupement par 
bandes : les composantes connexes de l'attracteur fusionnent successivement. 
On parle alors de cascade inverse. 

La complexité du comportement du système lorsque c varie provient de 
l'alternance de régimes périodiques et chaotiques. En effet, la situation pré- 
sentée avec une période 1 initiale se reproduit pour toutes les valeurs de 
périodes impaires. Ainsi, des fenêtres de périodicité s'installent brutalement 
à la suite de régimes chaotiques, s'achèvent par une cascade sous-harmonique 
suivie d'une cascade inverse et d'un régime chaotique. Ces fenêtres de pério- 
dicité sont denses dans [—2; 1/4], ce qui montre bien toute la complexité 
de la dynamique dans cette région de l'espace des paramètres. De plus, le 
complémentaire de cet ensemble ayant une mesure non-nulle, il reste pos- 
sible d'observer un comportement chaotique en choisissant le paramètre c 
aléatoirement suivant la mesure de Lebesgue. 

Une étude plus détaillée (et plus expérimentale) est faite dans |Man99] . 
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C.8 Comportement statistique des orbites 

L'étude des systèmes dynamiques mesurables est l'objet de la théorie 
ergodique, qui est notamment introduite dans [BP03| . Nous ne donnons ici 
que quelques définitions utiles pour notre étude. 

Définition C.31 (Mesure invariante). Une mesure fi est invariante par l'ap- 
plication (mesurable) / si pour toute partie mesurable A 

Définition C.32 (Ergodicité). Une application / : (X, i3, yu) — )■ {X,B,fï) 
qui préserve n est ergodique si 

\/AeB, r\A) = A^fiiA)e{0,l}. 

Proposition C.17. Une application f est ergodique si et seulement si toute 
application : X — t- C mesurable, telle que (p o f = (p presque partout, est 
presque partout constante. 

On peut remplacer dans cette proposition «mesurable» par L^{X,fi) ou 
L\X,^^). 

Exemple C.3. Les rotations d'angle a irrationnel et le doublement de l'angle 
sont ergodiques sur le cercle S^, pour la mesure de Lebesgue sur le cercle. 

La proposition suivante fait le lien avec la notion physique d'ergodicité. 

Proposition C.18. Si X est un espace métrique séparable, n une probabilité 
borélienne sur X , f : X X continue préservant X . Si f est ergodique, alors 
li-presque toute orbite est dense dans X . 

On a alors une estimation quantitative de la «densité» des orbites : pour 
toute partie mesurable A, la proportion de temps passée dans A par presque 
toutes les orbites est égale à fi{A). 

Théorème C.19 (Théorème ergodique de Birkhoff). Soit {X,B,fi) un es- 
pace mesuré, f : X —> X mesurable préservant fi. Pour tout (f) dans L^{X, fi), 
on note 

^ n—l 

Sn4>{x) = — 0(/'^(a;)) (somme de Birkhoff de (p). 

^ A:=0 

La limite 4>{x) = lim^^oo Sn4>{x) existe pour ^-presque tout x, cp o f = cp 
presque partout. Pour toute partie f -invariante A mesurable, de mesure finie, 
on a 

/ (pdjj, = / (pdfx. 
Ja J a 
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En particulier, si n est une mesure de probabilité ergodique, alors 




pour ^-presque tout x. 

Définition C.33 (Mesure physique, Ruelle-Bowen) . C'est une mesure de 
probabilité fi invariante par /, telle que pour toute application continue 
sur X, pour A-presque tout x G X, (f){x) = (pdfi. 

Cette condition est bien plus forte que l'ergodicité, puisque contrairement 
au résultat du théorème de Birkhoff, le résultat de convergence est vrai A-p.p. 
(A est la mesure de Lebesgue sur X, dont le support est X tout entier), et 
non /x-p.p., fi pouvant avoir un support bien moins grand que X tout entier. 
En particulier, si le support de /i a une mesure de Lebesgue nulle, le théorème 
ergodique de Birkhoff énonce un résultat que l'on n'observera jamais (p. s.) si 
l'on choisit une condition initiale x aléatoirement suivant A. De plus, lorsque 
la mesure physique existe (cela a été prouvé dans le cas du solénoïde), elle 
est unique (ce n'est pas toujours le cas pour les mesures ergodiques). 

La mesure physique (lorsqu'elle existe) donnant la densité de A-presque 
toute orbite, c'est elle que l'on observe empiriquement au cours des simula- 
tions numériques. 



C.9 Dimension fractale 

Certains des attracteurs que nous avons évoqués ont — au moins partiel- 
lement — une structure d'ensemble de CantorQj, de dimension non-entière. 
Nous allons donner un sens à cette affirmation, en définissant la dimension 
fractale d'un compact K. Il existe plusieurs autres notions de dimension non- 
entière (reliées les unes aux autres) , par exemple la dimension de HausdorffF^ ; 
la dimension fractale possède l'avantage d'être la plus simple à évaluer nu- 
mériquement. 

Définition C.34 (Dimension fractale). Soit K un compact d'un espace mé- 
trique {X,d). Pour tout e > 0, on note N^{K) le nombre minimal de boules 

46. Notamment le solénoïde, section [C.3.6[ et l'attracteur de Hénon, section [C. 3. 71 

47. on montre en général que la dimension de Hausdorfï HD{K) est inférieure où égale 
à Df{K). Il y a égalité pour des classes assez générales d'ensemble, par exemple pour 
l'exemple d'ensemble de Cantor décrit dans ce paragraphe. On pourra se référer à [PT93| 
pour le cas des ensembles de Cantor définis dynamiquement. 
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de rayon e nécessaires pour recouvrir K. La capacité limite ou dimension 
fractale de K est définie par 



Df[K) = limsup 

e^o - log e 

Lorsque K est une sous-variété de dimension finie, la dimension fractale 
est égale à la dimension topologique. 

Un autre cas classique est celui des ensembles de Cantor. Considérons un 
exemple où la dimension se calcule facilement : pour / = [a, b] intervalle, on 
note /(/) = [a,a + ^]U[a + b]. Le compact K = f]^^^ /"([0, 1]) est 

un ensemble de Cantor. On a N^n.[K) = 2", et e h-)- N^{K) est croissante, 
donc si 3" < e < 3"+\ 

nlog2 ^ log N,{K) ^ (n + l)log2 



(n + l)log3 —loge nlog3 
On en déduit que 

' log3 



C.IO Théorème de Whitney 

Pour visualiser les résultats des simulations numériques, nous avons pro- 
jeté en dimension 3 les points de MJ^ {N grand) que nous avions calculés. 
L'une des justifications a posteriori de la validité de la méthode est théo- 
rique et passe par le théorème de Whitney. En effet, l'attracteur semblant 
avoir une dimension (fractale) strictement inférieure à 1,5, il est possible de 
le plonger dans M^. 

Théorème C.20 (Whitney). Toute variété compacte lisse de dimensionn G 
N se plonge dans R^"^-'^. 

Ce résultat est démontré dans |Laf96] . Il se généralise au cas d'un compact 
de dimension fractale qui se plonge dans dès que > 2d. 



D Résultats détaillés 
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Figure 66 - Zone 1 : carte et points 01 à 04. 
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A =0.15 ipoint 5 [53^ 5^39-40] 



A =0.15 ;point 6 [62^ 6^39-40] 
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= 0.1 5 ; point 7 [ 62-64 -> 7 ^ 39-40 ] 



= 0.15 ; point 8 [ 52 ^ 8 ^ 30 ] 
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A|, = 0.15 ipoint 9 [51-75^9 ^40] 



Aj, = 0.1 5 ; point 10 [ 59-60 -f 10 ^ 68 ] 
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Figure 67 - Zone 1 : points 05 à 10. 
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Figure 68 - Zone 2 : carte et points 11 à 14. 
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A = 0.1 5 ; point 15 [ 39 ^ 1 5 ^ 28-1 6 



A =0.15 [point 16 [32-23^ 16 ^28-27 ) 




Ag = 0.15 ; point 19 [29^ 19 -s 46-34] 



Aj, = 0.1 5 ; point 20 [ 40 ^ 20 -s 26-27 ] 




Figure 69 - Zone 2 : points 15 à 20. 
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Figure 70 - Zone 3 : carte et points 21 à 24. 
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0.15 [point 25 [39-38^25 ^20-19 ) 



= 0.1 5 ; point 26 [ 12 ^ 26 ^ 13-12 ] 






0.15 : point27[ 20-20 -> 27^ 13-13 ) 



A =0.15 : point 28 [ 15-12 ^ 28 ^ 14-15 ) 




= 0.1 5 ; point 29 [ 07-06 -> 29 ^ 19 ] 



Aj, = 0.15 ; point 30 [ 08 ^ 30 ^ 18-19 ] 
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Figure 71 - Zone 3 : points 25 à 30. 
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Figure 72 - Zone 4 : carte et points 31 à 34. 
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A = 0.1 5 ; point 45 [ 17-18 ^ 45 ^ 42 



= 0.1 5 ; point 46 [ 19-18 ^ 46 ^ 21 ] 




Figure 75 - Zone 5 : points 45 à 50. 
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Figure 76 - Zone 6 : carte et points 51 à 54. 
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= 0.15 ; point 55 [41-42 (1^*) ^ 55 ^ 04-02 ] \ = « '^ : point 56 [ 67-80 ^ 56 ^ 08-03 ) 
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Ag = 0.15 ; point 59 [ 79 ^ 59 -s 10-09] Aj, = 0.15 ; point 60 [ 78 ^ 60 -s 10-72 ] 
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Figure 77 - Zone 6 : points 55 à 60. 
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Figure 78 - 



Zone 7 : carte et points 61 à 64. 
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A = 0.1 5 ; point 65 [ 77-78 ^ 65 ^ 72 ] 



= 0.1 5 ; point 66 [ 77-76 ^ 66 ^ 74 ] 




Figure 79 - Zone 7 : points 65 à 70. 
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Figure 80 - Zone 8 : carte et points 71 à 74. 
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A =0.15 ; point 75 [ 64 ^ 75 -f 76 ] 



= 0.1 5 ; point 76 [ 75 ^ 76 ^ 47-67 ] 




A = 0.1 5 ; point 77 [ 74 ^ 77 ^ 65-66 ] 



A =0.15 ; point 78 [ 73 ^ 78 ^ 60 ] 
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Figure 81 - Zone 8 : points 75 à 80. 
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Figure 82 - Localisation du pli et de sa préimage (deuxième méthode). 
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Courbure du pli {t=0) ; Maximum pourj=498 



Courbure du pli {1=1 .5) ; Maximum pour j=4E 





200 400 600 800 1000 1200 200 400 600 800 



(a) t = 

Courbur9 du pli (t=2) ; Maximum pour j=890 



(b) t = 1,5 

Courbure du pli (t=2.5} ; Maximum pour j=498 
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(c) t = 2 

Courbure du pli (1=3.3) ; Maximum pour j=890 




(d) t = 2,5 



Courbure du pli (t=4) ; Maximum pour j=8J 
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(e) t = 3,3 (f) i = 4 

Figure 83 - Formation du pli : courbure (deuxième méthode). 
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Formation du pli sur 4 ans 
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(b) Localisation du maximum 
Figure 84 - Formation du pli : maximum de courbure (deuxième méthode). 
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